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Abstract 

" In  order  to  investigate  the  heat  and  mass,  transfer  aspects  of 
therroomechanical  ablation,  ablation  experiments  were  performed  with 
models  made  of  pure  carbon  dioxide  (C0-)  and  a glass  bead -CO.,  composite. 

The  experiments  were  performed  in  a low  subsonic,  low  turbulence  j.et  at 

/ 

ambient  pressure  and  near  0?C  temperature.  The  model  geometry  was  a 
hemispherical  forebody-cylindrical  afterbody  with  a diameter  of  25  mn. 
The  stagnation  point  mass  transfer  of  the  pure  models  was  6.21  less 
than  predicted  by  equilibrium  theory.  The  distribution  around  the 
hemisphere  generally  agreed  with  the  predicted  laminar  heat  transfer 
result  of  monotonically  decreasing  heat  transfer  with  distance  from  the 
stagnation  point.  The  internal  temperature  response  was  approximately 
5%  higher  than  the  theoretically  predicted  wall  temperature.  These 
results  could  be  correlated  by  the  inclusion  of  a nonequilibrium 
.sublimation  model  into  the  theoretical  prediction  equations.  -The 
required  vaporization  coefficient  was  0.0005.  For  the  composite  models, 
the  stagnation  point  C02  mass  flux  was  23?  below  the  rate  of  the  pure 
CC>2  models.  The  mass  flux  around  the-  hemisphere  did  not  decrease  with 
distance  from  the  stagnation  point  until  beyond  the  45  degree  position, 
and  the  internal  temperature  was  slightly  higher  than  the  pure  CO^ 
models.  The  composite  results  could  be  predicted  by  reducing  the 
transfer  coefficient  by  6?  and  including  in  the  surface  model  diffusion 
of  the  sublimed  CO^  through  the  wells  between  the  beads,  nonequilibrium 
sublimation,  and  the  consequences  of  heat  transfer  through  the  beads  to 
the  subliming  surface. 
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THERMOMECHANICAL  ABLATION  ’ 

I.  Introduction 

The  reentry  of  a vehicle  into  the  atmosphere  cannot  be 
experimentally  duplicated  on  the  ground,  nor  can  the  vehicle 
in  flight  be  instrumented  well  enough  to  determine  all  of 
the  significant  phenomena  associated  with  nosetip  ablation. 
This  handicaps  the  analyst's  attempt  to  model  an  event 
governed  by  aerodynamics,  heat  transfer,  and  material 
science.  It  is,  therefore,  no  surprise  that  ablation 
modeling  is  a difficult  task  and  is  often  the  subject  of 
much  dispute.  . • 

Through  the  efforts  of  many  early  researchers,  the  rate 
of  thermochemical  ablation  can  be  approximately  determined 
with  rather  simple  relationships  which  yield  surprisingly 
accurate  predictions  for  some  materials.  The  key  to  the 
simple  solution  is  the  application  of  the  "Unity  Lewis 
Number”  assumption,  which  implies  equal  mass  and  thermal 
diffusivities,  to  the  surface  energy  and  specie  conserva- 
tion equations.  This  assumption  was  applied  to  the  problem 
by  Shvab  and  Zeldovich  (Ref  21  and  26)  m the  19^0 's  and 
received  its  widest  dissemination  when  presented  in  its 
most  general  form  by  Lester  Lees  at  the  Third  AGARD 
Colloquium  on  Combustion  and  Propulsion  in  1958  (Ref  12). 
Lees'  results  can  be  found  in  textbooks  such  as  Dorrance 
(Ref  3)  and  are  the  basis  for  most  of  the  ablation 


prediction  procedures  used  in  the  United  States  aerospace 
industry  today.  Spalding  (Ref  24)  employed  the  unity 
assumption  in  a slightly  different  presentation  to  yield 
equally  useful  results,  which  received  excellent  tutorial 
treatment  in  W.  M.  Kays’  text  (Ref  8).  These  formulations 
allow  the  determination  of  surface  temperature  and  the 
ratio  of  mass  transfer  rate  to  heat  transfer  coefficient 
when  the  boundary  layer  edge  conditions  and  the  thermo- 
chemistry of  the  system  are  specified.  Convective  heat 
transfer  theory  supplies  the  transfer  coefficient  which 
must  be  corrected  for  mass  addition  at  the  wall  (blowing). 
Kubota  (Ref  11)  simplified  lees’  theory  for  a binary  mixture 
with  nc  chemical  reactions  and  performed  some  early  low 

“ft 

temperature  ablation  experiments  using  ice  at  Mach  5*8. 

When  experimental  data  became  available  for  the 
ablation  of  graphite,  it  was  noted  that  the  mass  loss  rat';S 
were  higher  than  the  amount  predicted  by  thermochemical 
theory.  The  hypothesis  which  best  explained  the  data  was 
that  the  graphite  did  nrt  completely  transition  to  the 
vapor  phase,  but  that  solid  particles  were  removed  from  tlie 
surface  in  a high  enthalpy,  high  shear  ablation  environment 
(Ref  10).  In  fact,  Lundell  and  Dickey  were  able  to  photo- 
graph a shower  of  glowing  particles  downstream  of  an 
ablating  graphite  model  (Ref  13).  The  phenomena  is  called 
thermomechanical  ablation  or  erosion  and  it  was  estimated 
that  in  reentry  the  ratio  of  thermomechanical  ablation  to 
thermoehemical  ablation  for  graphite  could  exceed  one 


2 


(Ref  10) . As  the  understanding  of  turbulent  rough  wail 
heat  transfer  increased,  the  estimated  mechanical  loss 
fraction  was  revised  downward;  however,  uncertainties  exist 
due  to  the  fact  that  it  has  not  been  possible  to  directly 
measure  the  fraction  of  solids  removed  from  an  ablating 
graphite  surface. 

Ablation  performance  prediction  models  include  a 
provision  for  thermomechanical  mass  loss.  The  conservation 
of  mass  and  energy  equations  for  the  ablating  surface  allow 
solid  as  well  as  vapor  to  convect  through  the  boundary. 
However,  it  is  assumed  that  neither  the  surface  nor  the 
boundary  layer  phenomena  are  affected  by  the  solid  parti- 
cles (Ref  9 & 1 6). 

% ' 

The  objectives  of  this  research  are  to  perform  an 
ablation  experiment  in  which  a controlled  amount  of  solid 
material  is  released  from  an  abating  surface,  and  to 
improve  existing  thermomechanical  ablation  models. 

This  dissertation  contains  a derivation  of  the  approx- 
imate ablation  theory  which  was  used  to  select  solid  carbon 
dioxide  (CO,,  or  dry  ice)  from  among  various  candidates  for 
the  experiment.  The  construction  and  testing  of  both  pure 
COg  and  glass  bead-COg  models  are  described.  Following  a 
presentation  of  experimer: lal  results,  improvements  in  the 
thermomechanical  theory,  including  the  effects  of 
nonequilibrium  sublimation,  diffusion  wells,  and  microscale 
heat  transfer  events,  are  developed  and  discussed. 


3 


II.  Simplified  Ablation  Model 


The  primary  purpose  of  this  chapter  is  to  present  the 
approximate  ablation  theory  as  applied  to  a simplified 
model  of  a subliming  surface  including  theraomechanical 
ablation.  The  othe'r  analytical  relations  needed  to  deter- 
mine ablation  rates  are  also  presented.  A prediction  cf 
the  simplified  theory  is  then  compared  to  that  of  a more 
exact  solution. 

The  heat  and  mass  transfer  fluxes  are  shown  schemati- 
cally in  Figure  1.  The  coordinate  system  is  fixed  to  the 
solid-gas  interface.  The  system  is  assumed  to  have  only 
two  components:  the  subliming  substance  and  the  surround- 
ing gas,  "air.''  The  system  is  analyzed  at  a fixed  location 
of  the  streamwise  coordinate. 

The  overall  mass  flux  of  the  system  is  described  by 
the  equation 

(pv)w  + m*  = ip  (1) 

♦ 

where  nVj,  is  the  total  mass  flux  coming  to  the  surface,  m* 
is  the  mass  flux  which  leaves  the  surface  in,  the  solid 
phase,  and  (pv)  is ‘the  normal  mass  flux  of  the  gas  phase. 
If  DVj,  is  partitioned  into  the  amount  which  leaves  the 
surface  as  a solid,  m*,  and  the  amount  which  vaporizes,  m, 
the  overall  mass  conservation  equation  reduces  to 


Figure  1.  Thermochemical  - Thermoroechanical  Ablation  Process 


! 
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(2) 


(pv)w  = r,  ' 

The  mass  flux  of  the  vaporized  material  is  usually 
partitioned  into  convected  and  diffused  portions.  The 
convected  portion  is  expressed  as  the.  mass  fraction  of  the 
subliming  substance  at  the  wall,  K , times  the  normal  mass 

W 

flux  of  the  gas  phase  at  the  wall..  The  expression  for  the 
diffused  portion  is  derived  from  Pick's  law — the  diffused 
flux  is  proportional  to  the  concentration  gradient. 
Therefore,  the  conservation  of  mass  equation  for  the 
material  which  vaporizes  is 

(pv)wKw  - tpD  = m (3) 

The  equation  is  further  simplified  by  the  substitution 

of  (pv)  from  Equation  (2)  into  Equation  (3)  and  by 
w 

rearrangement . 

* 

i (1  - Kw)  = •-  (pD  §)„  CM 

The  energy  flux  of  the  system  is  described  by  the 
equation 

^rad^net  “ ^ 9y^w  “ ^-^PDhi  jfy^w  + ^v^whw 

+ **hs  + iint  = Vs  (5) 
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where  the  first  term  is  the  net  flux  of  heat  transferred 

by  radiation,  .the  second  term  is  the  classical  conductive 

term,  expressed  by  Fourier’s  Law  applied  to  the  gas  phase, 

the  third  term  is  the  enthalpy  carried  by  the  diffused 

portion  of  the  gas  phase  mass  flux,  and  the  fourth  term  is 

the  amount  carried  by  the  conveeted  portion.  The  amount  of 

mass  which  leaves  the  surface  in  the  solid  phase  is  assumed 

to  have  the  same  enthalpy  as  that  of  the  solid  coming  up  to 

the  surface,  h_  (Ref  10  & 1 6).  The  fifth  term  is  the 
s 

amount  of  heat  conducted  into  the  interior.  Subtracting 
m*hg  from  each  side  of  the  equation,  the  energy  equation 
becomes 


^rad-^net  " ^ 3y^w  '*  ^5^  ^pDhi  dy  \ 


•• 

+ <Pv)w>V  + = ”hS 


(6) 


Thus,  by  partitioning  the  total  mass  flux 'into  the 
amount  which  is  vaporized  and  that  which  is  not,  the 
equations  for  a system  with  thermomechanical  ablation  are 
transformed  into  the  same  equations  which  describe  a system 
with  pure  thermochemical  ablation:  Equations  (2),  (J), 

(*}•),  and  (6).  This  is  an  important  result  in  that  the 
development  which  follows  applies  to  either  pure  thermo- 
chemical ablation  or  combined  thermochemical  - thermo- 
mechanical ablation  provided  that  the  term  m is  restricted 
to  that  portion  of  the  mass  flux  which  is  vaporized  at  the 
surface.  This  result  also  applies  to  the  more  complicated 
ablation  process  of  graphite  (Ref  16),  and  is  the  basis  of 
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the  thermomechanical  ablation  experiment  of  this  research. 
Simply  stated,  the  existing  heat  and  mass  transfer  theory 
associated  with  thermomechanical  ablation  predicts  that  the 
response  of  the  vaporizing  material  does  not  depend  on  the 
amount  of  thermomechanical  ablation  in  progress.. 

To  simplify  the  following  development,  the  radiative 
heat  flux  term  and  the  internal  heat  conduction  term  are 
assumed  to  be  negligible  relative  to  the  other  terms  and 
are  eliminated  from  the  analysis.  Further  simplification 
of  the  energy  equation  results  from  the  assumption  that  the 
lewis  number  is  one . The  lewis  number  is  defined  as  the 
ratio  of  the  mass  diffusivity  to  the  thermal  diffusivity. 

• Le  = k7pC^  ’ (?) 

The  quantity  is  often  close  to  one  for  many  mixtures  of 
interest.  When  the  lewis  Number  equals  unity,  the  heat 
conduction  and  the  summation  of  the  diffusion  enthalpy  flux 
terms  of  Equation  (6)  combine  into  the  single  term 

This  combination  is  derived  in  Chapter  V.  The 
energy  equation  thsn  reduces  to 


m(h  -h  ) = (£-  %) 


8y  w 


At  this  point  it  is  convenient  to  replace  the  term 
with  the  derivative  by  a transfer  coefficient  times  a 
driving  potential.  The  most  convenient  coefficient  for  the 
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heat  and’  mass  transfer  problem  is  the  .Stanton  number, 
defined  by  the  equation 


C 


(9) 


' where  hr  is  the  recovery  enthalpy.  Substitution  into 
Equation  (8)  yields 

■<W  = (10) 

The  similarity  of  Equations  (4)  and  (8)  suggests  that 

* » 

t » , 

a transfer  coefficient  approach  might  also  be  useful  in 

Equation  (4).  In  fact,  with  no  chemical  reactions,  and 

♦ . 

■ when  the  Prand’tl  and  lewis  numbers  equal  unity,  the 

% 

•• 

boundary  layer  solutions  to  the  concentration  and  energy 
equations  are  similar,  which  will  be  shown  in  Chapter  V. 
With  similarity  the  identical  transfer  coefficient  can  be 
used  for  mass  transfer  with  the  concentration  difference 
replacing  enthalpy  difference  as  the  driving  potential. 

The  Stanton  number  for  mass  transfer  is 


C 


(11) 


Substitution  of  from  Equation  (11)  into  Equation  (4), 

with  K = 0,  yields 

© . 


K 


m 1 -WKW  PeueCH 


(12) 
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Equation  (10)  is  still  somewhat  difficult  to  use’  since 
the  enthalpy  of  the  gas  at  the  wall  is  a function  of  concen- 
tration and  temperature.  For  the  two  species,  (l)  the  sub- 
liming substance  and  (2)  "air”,  the  enthalpy  at  the  wall  is 

\ = Vlw  + (1-V  h2w  ' • 

Substitution  of  h_  from  Equation  (13)  into  (10)  yields 

w 

*whlw  + '1-Kw)  h2w-hs]  - (14) 

Pe^H^r  * Kwhlw  " (1"V  h2w^ 

Using  Equation  (12),  the  terms  of  Equation  (14)  can  be 
regrouped  such  that  all  terms  relating  to  the  enthalpy  of 
the  subliming  material  are  on  the  left  and  all  the  terms 
relating  to  the  enxhalpy  of  the  air  are  on  the  right. 

“ <hlw-hs5  = Peue°H.(Vh2w)  (15) 


The  enthalpy  difference  on  the  left  is  simply  the  enthalpy 
of  sublimation,  Ah 

Finally,  the  nondimensional  mass  transfer  parameter  B’ 
is  defined  by  the  equation 


B* 


m 


PeaA 


(16) 


) 

l 
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Substitution  of  B*  into  Equations  (12)  and  (15)  and 
rearrangement  yields 
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B' 


(1?) 


A 


B’ 


Ah 


sub 


(18) 


Assuming  that  the  gas  at  the  wall  can  be  modeled  as  a 
mixture  of  ideal  gases,  Equation  (17)  can  be  expressed  in 
terms  of  the  partial  pressure  of  the  subliming  material. 


B’  = 


M1  plw  / P 
m2  1 - (Plw/P) 


(19) 


Where:  M = Molecular  weight 

p,  - Partial  pressure  of  sub- 
x . liming  specie  at  the  wall 

* 

P = Mixture  pressure 

If  the  vapor  of  the  subliming  specie  is  assumed  to  be 
in  equilibrium  with  the  solid  phase,  the  partial  pressure 
is  the  equilibrium  vapor  pressure  which  is  a function  of 
the  wall  temperature. 

With  this  system  of  equations,  a unique  value  of  T 

w 

and  B'  can  be  determined  when  hr,  P,  and  the  thermochemis- 
try of  the  system  are  specified.  The  solution  of  these 
equations  can  be  obtained  graphically  by  assuming  various 
values  of  wall  temperature  and  plotting  B'  from  Equations 
(18)  and  (19).  Ablation  response  plots  are  shown  in 
Figure  2. • The  equilibrium  vapor  pressure  of  a material  is 
often  correlated  in  the  form 


log10(p)  = (-0.218 5 A/T)  + B . (20) 


where  p is  the  vapor  pressure  in  Torr  an*d  T is  the  tempera- 
ture in  degrees  Kelvin  (e.g.  The  Handbook  of  Chemistry  and 
Physics) . The  Clapeyron  equation  relates  the  enthalpy  of 
evapora-t-ion-or-subii-ina-tion— to-the-equ-i-l-i'briura-vapor- 


pressure — temperature  relation.  Equation  (20)  is  conven- 
iently arranged  so  that  the  constant  A is  the  molal  heat 
of  sublimation  Equations  (18),  (19) » and  (20) 

can  be  solved  by  an  iterative  process. 

Once  B’  and  T have  been  determined,  two  corrections 
w 

must  be  applied  to  conventional  heat  transfer -results 
before  mass  flux  can  be  determined.  First,  the  Stanton 
number,  which  is  based  on  the  actual  heat  transferred  to 
the  ablating  body,  must  reflect  the.  reduction  in  heat 
transfer  due  to  blowing.  Gazley,  Gross,  and  Masson  (Ref  5) 
reported  Equation  ,(21 ) , a correlation  of  theoretical 
results  for  the  ratio  o’f  Stanton  number  without  mass 
addition  to  the  boundary  layer,  C^Q,  to  that  with  mass 
addition. 


Mg  1//3 

1 + 0.64  (j^)  B’ 


(21) 


The  equation  compares  well  with  experiment,  especially  for 
low  B'  values. 

Second,  the  results  of  the  simplified  heat  and  mass 
transfer  equations  must  be  ccrrected  for  the  assumption  of 
unit  values  for  the  Lewis  and  Prandtl  numbers.  Analysis 


and  experiment  have  shown  that  moderate  departures  from 
unity  are  adequately  accounted  for  by  adjustment  of  the 
transfer  coefficient.  Kays  {Ref  8)  suggests  the  following 
lewis  number  correction  to  the  heat  transfer  Stanton  number 
for  mass  transfer  applications 

(Cu)  Mass  Transfer  0/~ 

—2 = (Le)  2/3  (22) 

(Oh)  Heat  Transfer 

In  summary,  the  mass  loss  is  estimated  by  the  follow- 
ing procedure: 

1)  B'  and  T are  determined  upon  specification 

w 

of  the  environmental  conditions  and  the  subliming  material 
using  a plot  or  by  otherwise  solving  Equations  (18)  and  (19)* 

2)  Cjj0/Cpj  is  determined  from  Equation  (21). 

3)  °Ho  is  determined  through  conventional  heat 

transfer  techniques.  Since  T has  been  determined, 

w 

reference  temperature  methods  can  be  used  for  properties. 

*0  Finally 

® = ^pe%CHo^S/CHo^Le^2^B'  ^23^ 

: • 

This  simple  theory  provided  the  information  necessary 
to  design  a sublimation  experiment;  however  a verification 
of  its  accuracy  was  required  before  it  could  be  used  to 
evaluate  experimental  results.  The  more  exact  solution  was 
obtained  using  the  Boundary  Layer  Integral  Matrix  Procedure 
(BLIMP)  written  by  Bartlett  and  Kendall  (Ref  2).  The 
BLIMP  solution  is  described  in  Appendix  A,  It  includes 
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the  mass  transfer  prediction  for  the  entire  hemispherical 
'forebody,  and  takes  into  account  the  individual  air 
species,  and  unequal  diffusion  coefficients.  A comparison 
* of  the  simple  theory  and  BLIMP  for  an  axisymmetric  stagna- 
tion point  for  the  approximate  conditions  of  the  experiment 
performed  in  this  research  is  shown  in  Table  I.  In  this 
case,  the  unity  Lewis  number  assumption  has  not  only 
simplified  the  derivation,  but  also  provided  very  accurate 
predictions  for  the  ablation  response.  The  boundary  layer' 
is  assumed  to  be  laminar. 

Table  I 

Comparison  of  BLIMP  and  Simple  Theory  Predictions  for  an 
Axisymmetric  Stagnation  Point 


BLIMP 

APPROXIMATE 

DIFFERENCE 

Tw  (K) 

181 

172 

-5.0* 

B' 

0.193 

0.192 

-0.5* 

PeueCHo 

0.262 

0.268 

2.3* 

CHo//CH 

1.094 

1.107 

0.9* 

peueCH 

0.239 

0.242 

1.3* 

m (kg/m  *s) 

•••”*  0.0461 

0.0465 

0.9* 

* 

H 


14 


III.  Experiment 


An  experiment  v/as  performed  which  provided  a direct 
comparison  of  ablation  with  and  without  solid  mass  addition 
to  the  boundary  layer.  This  chapter  contains  the  rationale 
used  in  the  selection  of  a subliming  material  for  the 
experiment,  the  method  used  in  the  fabrication  of  the 
models,  a description  of  the  experimental  test  facility, 
the  experimental  procedures,  and  the  methods  used  in  the 
reduction  of  the  data.  ■ > , 

Subliming;  Material  and  Environmental  Conditions 

* . 

The  initial  problem  was  to>find  a good  combination  of 

• 

subliming  material  and  environmental  conditions  for  the 
experiment.  Kubota's  paper.  (Ref  11)  was  quite  helpful  in 
this  regard  since  it  contained  graphs  of  ablation  response 
as  a function  of  environmental  conditions  for  water  ice, 
dry  ice,  and  camphor  which  are  summarized  in  Figure  2.  A 
low  turbulence  subsonic  free  jet  was  chosen  for  the  experi- 
ment because  of  its  ease  of  operation  and  low  cost.  The 
free  jet  afforded  open  access  to  the  model,  which  facilita- 
ted data  acquisition.  The  conditions  in  the  jet  were 
ambient  pressure  and  near  0°C  temperature.  With  this 
environment,  dry  ice  (solid  COg)  was  chosen  because  it 
responded  with  a B'  of  nearly  0.2,  which  is  high  enough  to 
yield  sufficient  mass  transfer  for  experimental  accuracy 
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Figure  2.  Ablation  Response  of  Dry  Ice,  Ice,  and  Camphor  (Summarized  From  Ref  11) 


and  is  within  the  range  of  B - values  experienced  by 
graphite  during  reentry. 

Axisymmetrical  flow  over  a hemispherical  forebody- 
* cylindrical  afterbody  was  chosen  over  a two-dimensional 
flow  configuration  for  a more  accurate  determination  of  the 
shape  change  during  the  experiment.  The  potential  flow 
solution  for  the  -selected  shape  is  one  of  the  many  computed 
by  A.M.O.  Smith  and  Pierce  and  presented  in  Reference  22. 
Their  computed  velocity  distribution  merges  into  that  for  a 
sphere  near  the  stagnation  point  as  shown  in  Figure  3* 
Therefore,  the  stagnation  point  heat  transfer  coefficient 
for  this  geometry  is  identical  to  that  of  a sphere.  - 

The  experiment  was  performed  at  Mach  numbers  below  0.1 
which  further  simplified  the  prediction  of  the  ablation 
response.  The  pressure  did  not  significantly  vary  around 
the  body  and  the  recovery  air  temperature  was  essentially 
constant  and  equal  to  the  free  stream  air  temperature. 

With  constant  environmental  conditions,  the  simple  equilib- 
rium theory  predicts  a constant  wall  temperature  and  B' 
around  the  body^,  which  further  implies  a constant  ratio  of 
Cjj  to  <W  The  mass  loss  is  then  directly  propoi-tional  to 

the  heat  transfer  coefficient  around  the  body  as  predicted 

, * 

by  Equation  (23).  This  conclusion  was  confirmed  by  BLIMP 
in  that  the  predicted  variation  of  wall  temperature  around 
the  body  was  only  0.11K  (0.0 65$)  and'  the  B'  variation  was 
only  0.0015  (0.?8#h 

m » 

v 

•* 
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Figure  3.  Boundary  Layer  Edge  Velocity  Normalized  by  Free  Stream  Vel 
and  a Hemisphere-Cylinder 
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Ablation  Models 

The  25  mm  diameter  models  were  produced  in  molds  of 
two  designs,  which  are  shown  in  Figure  4.  The  first  was 
made  from  seamless  tubing  with  a welded  end  plug  which  had 
a hemispherical  surface  machined  into  it.  A second  mold 
was  constructed  of  stainless  steel  tubing  with  a hemispher- 
ical surface  machined  into  one  of  the  threaded  end  plugs. 
The  second  design  was  much  easier  to  construct.  It 
facilitated  the  extraction  of  the  frozen  model,  since  both 
ends  of  the  tube  could  be  removed  and  the  model  pushed  out. 
Both  molds  had  a single  entranced , port  for  the  introduction 
of  carbon  dioxide. 

Thetfmomechanical  ablation  was  simulated  by  dispersing 
glass  beads  in  the  CO^.  Since,  the  speculated  size  of  the 
solid  particles  removed  from  ablating  graphite  was  on  the 
order  of  the  stagnation  point  bouniary  layer  thickness,  the 
boundary  layer  thickness  of  this  experiment  dictated  that 
120  {im  diameter  beads  be  chosen.  The  beads  were  screened 
to  improve  size  uniformity.  The  bead  shape  was  fairly 
spherical  as  shown  in  the  photomicrograph  in  Figure  5* 

The  most  reliable  method  of  obtaining  a uniform 
dispersion  of  glass  beads  In  the  C02  was  to  fully  pack  the 
beads  into  the  mold  before  introducing  the  C02.  -All  of 
the  glass  bead  C02  ablation  experiments  were  performed 
using  the  fully  packed  mode  of  glass  bead  loading.  A 
plastic  rod  (also  shown  in  Figure  4)  was  frozen  into  the 
center  of  each  model  to  serve  as  a solid  attachment  point. 
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a) 


The  rod  had  a, hole  drilled  along  its  axis,  so  that  a ther- 
mocouple junction  could  "be  placed  near  the  center  of  the 
hemispherical  nose. 

The  glass  head  and  C02  composite  models  were  quite 
sensitive  to  moisture.  If  any  moisture  was  present,  the 
beads  tended  to  agglomerate  during  the  experiment,  forming 
a crust  on  the  model  surface.  The  following  filling 
procedure  was  used  to  eliminate  the  moisture.  The  mold  was 
filled  with  preheated  beads,  then  maintained  at  3&0  K in  an 
oven  for  at  least  one  hour.  The  hot  mold  was  then  attached 
to  a system  for  introducing  COg.  into  the  mold  which  is 
shown  in  Figure  6. 


PRESSURE  V 


C02  > 
SUPPLY 
CYLINDER, 


REGULATOR 


CONNECTOR 


VALVE 


VACUUM 

PUMP 


PRESSURE 
G 


Figure  6.  System  for  Producing  Ablation  Models 


■A  vacuum  pump  attached  to  the  fixture  maintained  a very  low 
pressure  in  the  mold  while  it  was  cooling.  When  cool,  the 
mold  was  placed  in  an  insulated  container  and  packed  with 
dry  ice.  The  system  was  then  pressurized  with  CO 2 "to 
approximately  14  atm.  The  C02  pressure  was  maintained  for 
several  hours,  typically  overnight,  then  the  regulator  was 
closed.  Within  five  hours  the  pressure  in  the  mold 
decreased  to  slightly  above  atmospheric  pressure,  as  the 
remaining  C02  vapor  solidified. 

Experimental  Test  Facility  and  Procedures 

A free  jet  built  for  an  Air  Force  Institute  of  Technol- 
ogy' aeroaccoustic  study  (Ref  14).  supplied  dry,  very  low 
turbulence  air  to  the  models.  The  jet  exit  was  102  mm  x 
7 6 mm.  The  model  size  and  placement  were  consistent  with 
aerodynamic  testing  standards  required  to  simulate  free 
flight  through  still  air.  The  jet  and  model  are  shown  in 
Figure  7* 

Measurements  were  made  of  tunnel  plenum  temperature 
and  pressure,  ambient  pressure,  and  the  model  internal 
temperature  and  shape  as  a function  of  time.  Plenum 
pressure  was  measured  using  a precision  manometer.  The 
output  of  the  thermocouple  inside  of  the  model  was  recorded 
by  a continuous  strip  chart  recorder.  Model  shape  as  a 
function  of  time  was  obtained  from  35  mm  photographs  of  the 
model  with  a stop  watch  and  reference  scale  in  the  field  of 
view.  Typical  photographs  of  pure  C02  and  composite  models 
are  shown  in  Figure  8.  The  models,  positioned  in  front  of 


a flat  "black  background,  were  front  lighted  by  a strobe 
flash  attached  to  the  camera.  The  camera  was  a Mamiya/ 
Sekor  1000  TL  with  a 270  mm  lens.  With  this  system  the 
* camera  could  be  positioned  about  two  meters  from  the  model 
to  eliminate  parallax. 

Air  was  introduced  into  the  tunnel  at  the  desired  flow 
rate  and  the  ambient  conditions  were  recorded.  When  tunnel 
temperature  arid  pressure  stabilized,  they  were  recorded. 
Then  the  pressure  was  vented  from  the  mold,  the  mold  was 
opened,  and  the  model  was  extracted.  The  model  was 
attached  to  the  support  system  while  positioned  out  of  the 
jet  and  the  thermocouple  wires  were  connected  to  the. 
recorder  system.  The  support  was  then  swung  into  the  jet 
and  locked  into  position  on  the  jet  centerline.  Pictures 
were  taken  at  approximately  ten  second  intervals  until  the 
surface  receded  to  the  embedded  thermocouple. 

Data  Reduction 

All  of  the  data  required  to  determine  surface  reces- 
sion rate  are  contained  on  the  3 5 mm  film.  Each  frame  was 
projected  onto  a light  table  by  a precision  film  reader  and 
the  distance  reference  and  surface  profile  were  traced.  A 
typical  set  of  recession  profiles  is  shown  in  Figure  9. 

The  complete  set  of  recession  profiles  is  contained  in 
Appendix  B. 

The  conversion  of  the  recession  profiles  to  local 
recession  rates  around  the  hemispherical  surface  was 
facilitated  by  the  use  of  an  HP  9820A  Calculator  with  an 
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■input  digitizer  and  an  output  plotter.  Of  primary  interest 
was  the  normal  recession  of  the  stagnation  region  while  the 
surface  was  still  essentially  hemispherical.  In  prepara- 
* tion  for  the  data  reduction,  the  recession  profile  trace 
was  placed  on  the  digitizer  and  the  center  of  curvature  of 
the  front  surface  was  determined  with  the  aid  of  an  overlay 
of ' concentric  circles.  A protractor  was  then  placed  over 
the  recession  trace  with  its  axis  on  the  center  of  curva- 
ture and  with  the  center  ray  (labeled  zero)  perpendicular 
to  the  nozzle  exit  plane  as  shown  in  Figure  9*  The  radii 
of  the  protractor  are  15  degrees  apart. 

After  entering  the  time  associated  with  each  profile 
trace,  the  recession  along  each  radius  was  digitized.  As 
each  point  was  digitized,  the  location  was  stored,  the 
location  and  time  for  that  point  were  plotted,  and  the 
point  was  numbered  on  the  plot.  The  calculator  referenced 
all  distances  to  the  first  point  digitized  and  all  times  to 
the  time  associated  with  the  first  profile.  After  each 
point  along  the  ray  was  digitized,  points  were  selected  to 
include  in  a least  squares  linear  curve  fit.  The  equation 
of  the  line  was  then  printed  and  the  line  was  plotted  on 
the  recession  trace  through  the  range  of  interest.  The 
recession  rate  is  'the  slope  of  the  curve.  Two  typical 
recession  plots  are  shown  in  Figure  10.  The  essentially 
constant  stagnation  point  recession 'rate  shown  in  the  top 
plot  is  a' result  of  the  compensating  factors  of  decreasing 
•model  diameter  and  increasing  nose  bluntness.  In  the 
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second  plot,  the  curve  indicates  an  initial  rate,  associated 
^ith  the  hemispherical  shape,  which  decreases  as  model  shape 
change  decreases  the  transfer  coefficient.  All  of  the 
recession  plots  are  contained  in  Appendix  B. 

The  carbon  dioxide  mass  loss  rate  is  related  to  the 
normal  recession  rate,  s,  by  the  equation  = pCQ  s. 


The  density  of  solid  C02  in  the  temperature  range  of 

interest  was  determined  by  Maass  to  be  1.60  Mg/rtP  (Ref  15)* 
For  the  glass  bead  C02  mixture,  the  density  of  each 
constituent  was  determined  by  measuring  the  increase  in 
weight  of  the  mold  after  addition  of  the  glass  beads,  then 
after  C02*  solidification.  The,  center  rod  was  not  used  in 
the  five  models  especially  fabricated  for  density  determi- 
nation. The  density  of  glass  in  the  mixture  was  1.48  ± 
0.02  Mg/iP  of  mixture  and  the  density  of  the  composite  was 
2.09  * 0.03  Mg/m^.  Therefore,  the  density  of  the  C02  was 
0.6l  - 0.05  Mg/rP  of  mixture. 


Experimental  Conditions  and  Physical  Properties 

The  range  of  experimental  conditions  and  physical 
properties  of  the  experimental  materials  are  listed  in 
Table  II.  The  glass  bead-C02  composite  thermal  properties 
are  based  on  the  properties  of  the  constituents.  Since  the 
conductivities  of  the  materials  are  nearly  identical,  the 
effective  conductivity  of  the  composite  is  essentially  that 
of  the  continuous  phase,  C02,  0.71  W/m*K  (Ref  7).  The 
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density  of  the  composite  was  determined  to  be  2.09  Mg/nP. 

By  adding  the  mass  weighted  contribution  of  the  specific 
heat  of  each  material,  the  specific  heat  of  the  composite 
is  determined  to  be  0.92  kJ/kg*K.  The  effective  thermal 
diffusivity  of  the  composite  is,  therefore,  0,37  mm^s. 

This  result  is  quite  fortuitous,  since,  with  nearly  identi- 
cal thermal  conductivity  and  diffusivity,  the  pure  CO^ 
models  and  the  composite  models  respond  essentially  the 
same  in  both  transient  and  steady  state  conduction. 
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Table  II  . 


Range  of  Experimental  Conditions 
and  Model  Component  Properties 


MODEL  RADIUS  (mm) 

10.0  - 

12.7 

VELOCITY  im/s) 

18.6  - 

38.3 

AIR  TEMPERATURE  (K) 

270  - 282 

PRESSURE  (atm) 

0.957  - 

- O.98O 

REYNOLDS  NUMBER  ReR 

19,500 

- 32,400 

PeueCH  Ocg/ro2*  s) 

0.21  - 

0.33 

PRANDTL  NUMBER  Pr* 

0.735 

LEWIS  NUMBER  Le* 

0.886 

* 

SCHMIDT  NUMBER  Sc* 

0.829 

SOLID  PROPERTIES  AT  170K 

CARBON 

DIOXIDE 

GLASS 

COMPOSITE 

DENSITY  (Mg/m3) 

1.60  (15)** 

2.7  (4) 

2.09 

THERMAL  CONDUCTIVITY 
(W/m-K) 

0.71  (6) 

0.79  (4) 

0.71 

SPECIFIC  HEAT  (kJ/kg*K) 

1.2  (15) 

0.84  (4) 

0.92 

THERMAL  DIFFUSIVITY 

(mm2/s) 

0.38 

0.35 

0.37 

MASS  FRACTION  IN 
COMPOSITE 

0.29 

0.71 

ENTHALPY  OF  SUBLIMATION 
(kJ/kg) 

589.1  (6) 

C02  VAPOR  PRESSURE  Log 

1QP  (atm)  = 

-1352.4/TOO 

+ 6.9576  + 

* At  Eckert's  reference  temperature 

(approximately  222K) 

**  Number  in  parenthesis  is  the 

reference  number 

+ This  correlation  is  based  on  the  enthalpy  of  sub- 
limation at  170  K and  the  vapor  pressure  data 
reported  in  Reference  6. 
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IV.  Experimental  Results 

The  mass  loss  results  for  the  pure  COg  models  are 
shown  in  Figure  11.  The  data  have  been  normalized  by  the 
stagnation  point  mass  flux  predicted  by  BLIMP,  adjusted  to 
the  particular  conditions  of  the  experiment  by  the  approxi- 
mate theory  presented  in  Chapter  II.  The  normalized  BLIMP 
distribution  is  also  shown.  Except  for  the  stagnation 
point  value,  each  data  point  shown  is  the  average  of  the 
values  on  either  side  of  the  stagnation  point.  The  results 
are  in  fair  agreement  with  the  predicted  resuits s although 
on  the  average,  the  experimental  values  are  less  than 
predicted,  -6. 5#  less  at  the  stagnation  point. 

The  normalized  carbon  dioxide  mass  transfer  results  of 
the  glass  bead  laden  models  are  shown  in  Figure  12.  The 
average  stagnation  point  carbon  dioxide  mass  loss  rate  was 
28#  less  than  the  amount  predicted  by  the  approximate 
theory,  or  23#  less  than  the  experimentally  observed  mass 
loss  rate  of  the  pure  carbon  dioxide  models.  The  distribu- 
tion of  mass  loss  around  the  hemisphere  was  different  from 
the  pure  C02  results  in  that  the  mass  loss  rate  did  not 
decrease  until  beyond  the  45  degree  location.  The  mass 
loss  data  for  these  models  also  exhibited  more  scatter, 
which  can  be  attributed,  at  least  in  part,  to  the  variation 
of  the  COg  mass  fraction  in  the  composite  models. 
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Figure  11.  CO2  Mass  Loss— Experimental  Data  From  Pure  CO2  Models 

and  Theoretical  Prediction— All  Normalized  by  Theoretical 
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The  thermal  response,  as  sensed  by  the  embedded 
thermocouple,  of  all  models  can  be  divided  into  three 
periods.  Referring  to  the  labels  on  the  example  tempera- 
ture trace  shown  in  Figure  13 , there  was  a period  of 
decreasing  temperature  (A  to  B),  a period  of  approximately 
constant  temperature,  referred  to  as  the  asymptotic  temper- 
ature (B  to  C),  and  a period  of  increasing  temperature  (C 
to  the  end  of  the  test) . The  initial  temperature  drop  is 
an  indication  that  the  temperature  of  the  dry  ice  packed 
container  in  which  the  mold  was  placed  while  the  model  was 
solidifying  was  higher  than  the equilibrium  wall  tempera- 
ture during  the  convective  heat  and  mass  transfer  process. 
The  increase  in  temperature  indicates  a response  to  possi- 
ble  thermal  gradients  near  the. Surface  and,  upon  exposure, 
the  warmth  of  the  air  stream. 

In  a few  experiments,  the  thermocouple  data  was 
obviously  erroneous,  in  others,  the  mass  loss  data  was 
useless  because  of  flaws  in  the  model,  yet  the  temperature 
data  was  acceptable.  A table  of  asymptotic  temperatures  is 
included  in  Appendix  B.  Of  three  experiments  with  pure 
COg,  the  mean  asymptotic  value  was  179*3  K.  Of  ten 
experiments  with  composite  models,  the  mean  asymptotic 
value  was  182.1  K. 
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V.  Ablation  of  the  Composite  Material 

The  addition  of  nonsubliming  material  to  the  CO^  had  a 
significant  effect  on  the  observed  mass  transfer  from  the 
model,  A change  in  mass  transfer  could  be  attributed  to  a 
change  in  either  the  transfer  coefficient,  pnu  C„  or  the 
driving  potential,  3'.  Had  the  reduction  in  mass  loss  been 
solely  the  result  of  a reduced  transfer  coefficient,  no 
change  in  model  temperature  would  be  expected.  The 
slightly  higher  model  temperatures  associated  with’  the 
composite  material  and  the  high  percentage  of  nonsubliming 
material 4 dictated  that  a more . thorough  analysis  of  the 
ablating  surface  be  performed..  ‘ 

Inspection  of  the  enlarged  speculated  model  of  the 
surface  shown  in  Figure  14  suggests  items  for  consideration. 
First,  the  wall  will  most  certainly  be  rough  rather  than 
smooth.  On  a discrete  level,  the  raised  glass  bead  surface 
probably  receives  more  heat  flux  than  the  lower  C02  surface, 
especially  considering  the  fact  that  the  C02  is  outgassing, 
while  the  glass  is  not.  The  glass  beads  affect  the  mass 
transfer  of  the  C0?  by  constricting  the  surface  area 
available  for  sublimation,  which  causes  a higher,  local  mass 
transfer  rate  at  the  subliming  surface.  The  presence  of 
the  beads  also  caused  the  C02  vapor  to  diffuse  through 
wells  between  the  beads  to  the  boundary  layer. 
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Figure  14-.  Model  of  Ablating  Glass  Bead  - CO^  Surface 

In  this  chapter  the  heat  and  mass  transfer  equations 
which  were  presented  in  simplified  form  in  Chapter  II  are 
reexamined  with  special  attention  given  to  the'  probable 
microcosm  of  the  glass  bead-C02  surface.  The  improved 
■equations  are  then  individually  investigated,  first  the  new 
conservation  of  energy  equation,  then  the  conservation  of 
mass  equation.  Finally,  the  computational  scheme  used  to 
solve  the  coupled  equations  is  presented. 

Improved  Ablation  Model 

To  improve  the  modeling  of  the  ablating  surface,  the 
single  control  surface  is  replaced  by  the  control  volume 
shown  in  Figure  15.  The  top  of  the  control  volume  is  the 
fluid  dynamic  wall,  which  can  be  roughly  identified  as  the 
top  of  the  glass  bead  layer.  The  bottom  of  the  volume  is 
the  subliming  surface. 

The  total  solid  mass  flux  entering  the  control  volume 
is  immediately  partitioned  into  that  which  will  sublime,  m, 
and  that  which  will  not,  m*.  The  expansion  of  the  control 
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q = heat  conducted  into  the  interior 
int 

^rad^net  = net  radiative  heat  flux 


= mass  fraction  of  the  i species 
p = density  of  the  gas 
v = normal  velocity  of  the  gas 
h = .enthalpy 

m ~ mass  flux  of  the  subliming  material 

y = distance  above  surface 

T = temperature  of  the  gas 

D = binary  diffusion  coefficient 

k = thermal  conductivity  of  the  gas  mixture 

m*  = glass  bead  mass  flux 

• h*  = enthalpy  of  the  glass  beads  at  TQ 

Ah*  = increase  in  enthalpy  of  beads  through 
control  volume 
% 

Subscripts: 

s = solid  phase  immediately  below  subliming 

surface 
, • 

, *« 

w = gas  phase  at  fluid  dynamic  wall 

1 = subliming  specie 

2 = "air" 

Figure  15-  Control  Volume  for  Ablation  Model 
of  the  Composite  Material 
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surface  to  a control  volume  does  not  affect  the  conserva- 
tion of  overall  mass  and  individual  specie  equations. 
Rather  than  repeat  the  previous  development,  the  important 
results  are  listed. 


• 

m 

’ (Py)W 

(2) 

m(l-K„) 

” ff>w 

(4) 

The  improved  conserv  tion  of  energy  equation  contains 
all  of  the  terms  identified  in  Chapter  II,  including  the 
radiative  heat  flux  term  and  the  internal  heat  conduction 
term,  and  some  additional  terms  associated  with  the  flux  of 
the  glass  beads.  The  energy  terms  are  shown  schematically 
and  defined  in  Figure  15.  The  energy  equation  for  this 
system  is 


ihs+mV8  - -[<*“>, 

+ i*(h;  + w*)  + (4rad)net 


+ j,  <e6hi 

+ ( PV)W  + ‘lint 


(24) 


The  terms  in  the  bracket  may  be  simplified  under  the  unity 
Lewis  number  assumption.  Since  this  is  a binary  mixture, 
the  diffusion  coefficient,  D,  is  the  same  for  both  species 
(Ref  4s 4505.  First  the  Lewis  number,  which  was  defined  in 
Equation  {7),  is  substituted  into  the  bracketed  expression. 
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(25) 
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By  definition,  the  specific  heat  of  a nonreacting  mixture  is 


C = Z C . K. 
P i=1  Pi  i 


(26) 


and,  from  thermodynamic  analysis,  the  specific  heat  at 
constant  pressure  of  a specie  is 


3h. 

cPi  = hnPp 


(27) 


Substituting  these  definitions  into  Equation  (25),  the 
expression  becomes 

/ k \ r — ah*  <\tn  2 3K* 

(cp  w Ci=i K.  jT  + hi  sA  (28) 

Given  the  general  case,  that  enthalpy  is  a function  of 
temperature  and  pressure,  the  derivative  of  enthalpy  with 
distance  y is 


3h 

• ay 


r 3h\  9T  , / 3h\  aP 
' 3T ' p 3y  v 3P ; T ay 


(29) 


The  enthalpy  of  an  ideal  gas  is  a function  of  temperature 

only,  but  for  any  gas  the  second  term  is  zero  since  the 

% • 

pressure  is  assumed’ to  be  constant  across  the  boundary 
layer  in  the  normal  direction.  Therefore  Equation  (28) 
Simplifies  to 


k _2  ah*  2 gK* 

%>„  Ki  3F  + I-  if1  ^ 3^]w  (30) 
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Now,  if  it  is  assumed  that  Le  = 1,  the  sum  can  be  expressed 
as 


v 2 9h-  3K.  ' 

^Ki  "ay”  + hi  3y 


(31) 


Combining  the  differentials,  the  term  simplifies  to 


{C \ 3y  Kihi^w 


(k_  ih) 

vCp  3y'w 


(32) 


The  combined  term  is  often  referred  to  as  the  heat 
conduction  term  and,  as  defined  5.n  Chapter  II,  it  is 
expressed  as  a coefficient  times  an  enthalpy  difference 

l^w  = peueCH  (W  (33) 

x* 


The  driving  potential  is  the -difference  between  the 
recovery  enthalpy  and  the  enthalpy  of  the  mixture  at  the 
wall.  The  Stanton  number  is  substituted  into  Equation  (24), 
and,  after  rearrangement,  the  equation  becomes 


peueCH  (hr-hw}  " " ™*Ah* 

^rad^net  + ^-int 


(34) 


Further  simplification  requires  a boundary  layer 
similarity  result.  The  steady  state  conservation  of  mass 
equation  for  a specie  in  a nonreacting  boundary  layer  is 
(Ref  12) 
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(35) 


where s s = the  streamwise  coordinate 
u = streamwise  velocity 


For  the  same  system,  the  total  enthalpy  form  of  the 
energy  equation  is 


3H  4.  9H  - 3 / vi  3H  , 9 r m 1 ) a(u 
pu35+pv3^  - ^ sP  ay~ 


* ^ 3K’- 

+ 4CpD!l-feV?.  hi3F^ 


(36) 


where:  H = h + 


Assuming  that  the  Prandtl  'and  lewis  numbers  both  equal 
one,  the  last  two  terms  of  the  energy  equation  drop  out. 
This  assumption  implies  that  the  Schmidt  number  (Sc  = pD/y), 
the  ratio  of  Prandtl  number  to  Lewis  number,  is  also  one. 
Therefore,  m can  be  substituted  for  pD  in  Equation  (35)* 
and  the  equations  become 


3K.  3K. 

PU  3ii+  PT5T 


__  d / 1 \ 

3y  3y  5 


(37) 
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One  solution  for  Equation  (37)  and  (38)  is 


Kl  = AjH  + Bj_ 


(39) 
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where  A^  and  B^  are  constants  for  a given  specie.  The 
constants  for  the  subliming  specie  are  evaluated  by 
applying  -the  following  boundary  conditions: 


At  the  fluid  dynamic  wall: 
= Kw  and  H = H„  = h„  + 


At  the  boundary  layer  edge: 


K = K = 0 
e 


H = H 


Equation  (39)  then  becomes 


. Ky, 

■ (VV 


(40) 


The  derivative  of  Equation  (40)  with  respect  to  y, 
evaluated  at  the  wall  is  ,, 


r3K>  _ 

(Vw 


e w 


(®) 

v 3y'w 


(41) 


and  since  H = h + u /2,  then 


+ t/ 

' 3y  w 3y  w / w 3y  i 


(42) 


Substitution  of  Equation  (42)  into  (4l)  yields 


(3K) 

' 3y'w 


- --w  - (J&) 
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(43) 


When  Equation  (43)  is  substituted  into  the  conserva- 
tion of  specie  equation,  Equation  (4) , the  result  is 


m 


PD  Kw  /Jh% 
n-Kw)  (He-hw)  1 3y  w 


(44)- 


The  total  enthalpy  at  the  edge,  H , is  equal  to  the 

v 

recovery  enthalpy,  h , either  at  the  stagnation  point,  or 
everywhere  if  the  Prandtl  number  is  one.  Application  of 
the  unity  Lewis  number  assumption  and  the  definition  of 
the  Stanton  number  finally  yields  Equation  (12) 

m = ^1-K~^.  peueCH  {12) 

W 

As  in  Chapter  II,  the  energy  equation  -is  now  re- 
arranged using  the  definition  of  the  enthalpy  of  a mixture 
and  Equation  (12)  to  yield 


peue¥] hr'h2w)  ~ ' **Ah*  = 

^rad^net  + ^int 


(45) 


Since  the  fluid  dynamic  wall  temperature,  T , can  be 

w 

different  from  the  temperature  at  the  subliming  surface, 

T , the  enthalpy  difference  (h,  - h_)  can  be  divided  into 

8 ^ J.W  S 

the  enthalpy  of  sublimation  at  T , A h plus  the  enthalpy 
required  to  raise  the  temperature  of  the  vapor  to  Tw>Ahw  g. 
When  this-  distinction  and  the  definition  of  B ' are  applied 
to  Equation  (45),  the  result  is 
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fv  <vV  - i(ihSHh  + Ahw>s)--  ;*4h*  = 

^rad^net  + ^-int 


(46) 


Solving  for  B',  Equation  (46)  becomes 


B’  = 


(hr-h2w) 


<ihsub  + Ahw,s  + £ ^ + + — ) 


.(47) 
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If  each  term  in  the  denominator  is  normalized  by  Ahgut),  the 
result  is 

h - h. 


B'  = 


*2w 


"Vb'1  + x> 


(48) 


where:.  X = E X. 

i=l  x 


and 


xi  = 


rad  net 


mAh 


sub 


7 = m*&h^  tr  _ A*V.s  y _ ^int 
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The  nondimens ional  parameter  X is  the  sum  of  the  relative 
contribution  of  each  term  which  was  neglected  in  Chapter  II. 

The  definition  of  B'  is  again  applied  to  the  continu- 
ity of  specie  equation  which  can  either  be  expressed  in 
terms  of  mass  fraction  or  partial  pressure. 


B' 


Kw 

1 -Kw 


(1?) 


M1  Plv/P 
= m2  1 - lPlw/PT 


(19) 


1*6 


While  Equations  (17)  and  (19)  are  identical  to  those 

derived  in  Chapter  II,  neither  the  concentration  nor  the 
• • 
partial  pressure  will  he  assumed  to'  he  equilibrium  values. 

The  factors  which  influence  the  evaluation  of  the  concen- 
tration at  the  wall  will  he  discussed'  following  the 
modeling  of  the  additional  terms  of  the  energy  equation. 


Equation 


The  nondimensional  parameter  X in  the  energy  equation 
can  have  either  a positive  or  a negative  value.  If  X were 
zero,  Equation  (18)  of  Chapter  II  would  he  generated.  Each 
contribution  to  X will  he  separately  evaluated. 


Radiation  Heat  Transfer.  X^  is  the  ratio  of  net  heat 
transfer  by  radiative  process  to  the  rate  of  heat  absorp- 
tion due  to  sublimation.  Assuming  the  surroundings  are 
radiating  as  a black  body  at  room  temperature,  the  net  flux 
can  be  approximated  by 


...  II 

(q  ,)  , = EOT  - a of 

Hrad  net  w room 


(49) 


Further,  assuming  that  the  temperature  difference  is  small 
enough  that  e and  a are  approximately  equal,  the  equation 
simplifies  to 


iHrad  net 

— 1 

m Ah_  , 


/ 4 4 ft 0) 

<Tw-Troom)  15 


4? 


A pure  C02  surface  and  a glass  iaden  C02  surface  would 
probably  have  different  e's,  however  both  are  assumed  to  be 
unity  for  this  approximation. 

The  next  two  contributions  to  X reflect  the  fact  that 
some  energy  can  be  absorbed  by  the  solid  and  the  gas  pass- 
ing through  the  control  volume.  The  evaluation  of  these 
terms  requires  an  estimate  of  the  temperature  difference 
across  the  control  volume.  In  the  development  which 
follows,  a sphere  is  analyzed  to  determine  the  temperature 
rise  that  is  caused  by  an  imposed  heat  flux,  and  the 
increase  in  energy  implied  by  that  temperature  rise.  Next 
an  analysis  is  performed  to  determine  the  thermal  resist- 
ance of  the  gas  between  the  spheres.  An  analysis  of  the 
relative  resistances  of  the  sphere,  the  gas  columns,  and 
the  boundary  layer  finally  yields  an  estimate  of  the 
temperature  drop. 

Sphere  Analysis.  It  will  be  shown  that  the  majority 
of  the  heat  flux  to  the  model  passes  through  the  glass 
beads.  The  thermal  analysis  of  the  beads  requires  certain 
assumptions  which  cannot  be  directly  verified.  It  is 
assumed  that  the  beads  are  spherical,  and  that  they  do  not 
leave  the  surface  until  the  C02  surface  has  receded  to  the 
bead  equator.  It  is  further  assumed  that  the  thermal 
response  of  the  sphere  is  symmetrical  about  the  axis  normal 
to  the  model  surface  - temperature  is  a function  of  radial 
position  r,  and  angle  9,  Figure  16.  While  it  might  be 
suspected  that  transient  heat  conduction  solution  would 
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.be  required  to  obtain 
the  thermal  response 
of  the  sphere,  steady- 
state  solutions  are 
adequate.  Even  at  the 
highest  mass  transfer 
rate,  the  surfaces* 
receded  at  about  0.07 
mm/s,  traveling  one 
bead  depth  is  about  • 
two  seconds.  The 


Figure  16.  A Sphere  Half  Emerged 
From  the  COg 


Fourier  modulus,  defined  as  the  thermal  diffusivity  times 
the  time  divided  by  the  square  of  the  radius,  is  the 
nondimensional  time  scale  with  which  transient  fluxes  are 
correlated.  By  the  time  the  modulus  is  equal  to  unity,  the 
transient  flux  in  a sphere  is  negligible.  With  the  small 
radius  of  the  glass  beads,  the  modulus  becomes  greater  than 
one  within  0.01  seconds.  The  steady  state  heat  conduction 
equation,  in  spherical  coordinates,  for  the  axisymmetrical 
temperature  distribution  is  (Ref  20:135) 


2r  U ♦ r2  9'T 


3r 


3 r2 


+ ctn  8 U + ~^2 
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3 0' 


= 0 


(51) 


There  is  an  analytical’  solution  to  this  problem  in  the  form 
of  the  Legendre  polynomial  series.  Unfortunately,  unless 
the  surface  temperature  is  of  a fortuitous  distribution, 
convergence  can  require  a prohibitive  number  of  terms. 


While  several  numerical  programs  were  available  which 
could  solve  the  conduction  problem  in  rectangular  form,  it 
was  decided  that  it  would  be  more  fruitful  to  write  a 
relaxation  routine  to  solve  the  spherical  problem  directly. 
The  program  computed  internal  temperatures  for  a specified 
surface  temperature  distribution,  the  heat  flux  through  the 
sphere,  and  the  amount  of  energy  stored  in  the  sphere  above 
the  amount  it  would  have  had  at  a uniform  minimum  tempera- 
ture level. 

Solutions  were  obtained  for  two  cases:  I)  a cosine 

distribution  from  a maximum  T-yj  at  0 = 0 to  the  minimum 

T at  0 = it,  II)  a cosine  distribution  over  the  top 
o 

hemisphere  and  a constant  Tq  on  the  bottom  hemisphere.  The 
first  case  is  actually  the  first  Legendre  polynomial,  which 
supplied  a good  test  case  for  the  program.  The  resulting 
temperature  distributions  are  shown  in  Figure  17  and  the 
nondimensional  heat  flux  and  energy  storage  are  shown  in 
Table  III. 

These  two  cases  bracket  the  expected  result.  If  the 
temperature  distribution  is  symmetrical  about  the  equator, 
as  in  Case  I,  the  mean  energy  stored  in  the  bead  is  half 
the  amount  it  would  have  if  it  reached  a uniform  T^.  If 
the  C02  maintained  the  bottom  surface  at  an  approximately 
constant  temperature,  Tq,  the  energy  gain  is  about  one 
fourth  the  uniform  value.  The  fluxes  remained 
essentially  equal.  These  results  were  the  relations  used 

to  estimate  the  expected  temperature  difference  and  the 
resulting  increase  in  energy  because  of  that  difference. 
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Table  III 


Summary  of  Theoretical  Predictions  for  Heat  Flux 
and  Energy  Gain  of  a Sphere 

AE 

pSp  CV3)Ta3(TM-T0) 

0.5 
' 0.24 

Gas  Column  Analysis.  In  order  to  estimate  the  heat 
flux  through  the  bead,  the  thermal  resistance  of  a gas 
column  must  be  estimated.  For  analysis,  assume  the 
idealized  geometry  of  three  spheres  half  exposed  from  the 
surface  as  shown  in  Figure  18.  The  column  of  interest 
begins  as  a triangle  at  the  top  of  the  beads  with  an  area 
of  a and  bottoms  in  the  shaded  region  which  has  an 

2 ••  7T  2 

area  2~— — a • The  area  at  any  intermediate  plane  x is 

AM  = f C + x2]  (52) 

% 

Assuming  no  heat  flo w out  of  the  gas  column  to  the  glass 
walls,  a quasi-one  dimensional  expression  for  the  heat 
transferred  through/ahy  plane  is 

Q = - k ^ A(x)  (53) 


CASE 


klTa(TM-To) 


I 

II 


0.5 
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The  solution  of  this  equation  with  the . boundary  conditions 

that  the  temperature  at  the  top  of  the  column  is  T ^ and  at 

the  bottom  is  T is 
o 


2 Q 
kua 


(5*0 


Evaluating  the  constants  and  rearranging  the  equation 


3delds 


Q . 

* ka  ITm  - Tq) 


0.127 


(55) 
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Figure  18.  Idealized  Gas  Column  Geometry 


Resistance  Network  Analysis. . Finally  a -network  can  be 

constructed  to  represent  the  system.  A "unit  cell”  of 

surface  consists  of  one  glass  bead-  and  two  gas  columns  with 

__  2 

a total  cross  sectional  area  of  2 /J  a . The  network  is 
shown  in  Figure  19 • The  heat  flux  through  the  boundary 
layer  is  v, 

Q 


‘BL  = P.Wp'VV4 


Therefore,  the  thermal  resistance  is 


(56)  • : M 
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'EL  Q 


BL 
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(57) 


For  the  glass  bead,  using  the  results  of  the  sphere  analysis 


«GB  = ^ * *GB  *<.VV 


(58) 


Therefore 


GB  = 


^GB 


(59) 


For  two  gas  columns 


9o  “ 2«>^AS  a<Vts,(0-127) 


(60) 


Therefore 


R 


"^GAS 


(61) 


The  overall  resistance  of  the  network  is 


R = Rtjt  + 


RGB  RC 


BL  Rqb  + Rq 


5^ 


The  resistance  of  the  glass 
head  is  about  two  orders  of 
magnitude  less  than  either 
* the  boundary  layer -or  the 
air  columns;  therefore s the 
overall  resistance  to  heat 
flux  is  essentially  the 
boundary  layer  resistance. 

Since  the  temperature 
difference  imposed  on  the 
columns  and  the  bead  is  the 
same 
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Figure  19-  Resistance 
Network 
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The  ratio  of  the  heat  fluxes  is 


9 


Considering  the  order  of  magnitude  of  the  resistances, 
essentially  all  of  the  flux  is  through  the  beads.  The 
flux  through  the  bead  is,  then,  approximately  equal  to  the 
flux  through  the  boundary  layer;  therefore 


AT. 


BL 


R 


BL 


AT, 


GB 


R 


GB- 


(65) 


Equation  (65)  is  the  important  relation  for  estimation 
of  the  temperature  difference  between  the  fluid  dynamic 
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wall  and  the  subliming  surface,  which  is  assumed  to  be  • 

Enthalpy  Increase  of  the  Evolved  Solids.  The  enthalpy 
increase.  Ah  , corresponds  to  the  energy  increase  per  unit 
mass  as  predicted  by  the  glass  bead  analysis. 


A E 

P (4/3  * a?  ) 


(66) 


It  is  assumed  that  the  actual  value  of  this  quantity  is 
bracketed  by  the  predictions  for  Case  I and  Case  II  shown 


in  Table  III. 
then  X2  is 


If  the  average  of  the  two  cases  is  used, 
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Enthalpy  Increase  of  the  Evolved  Gases.  is  the 
ratio  of  the  heat  absorption  rate  .of  the  gas  while  increas- 
ing in  temperature  from  T to  T , to  the  heat  absorption 

s w 

rate  due  to  sublimation,  which  is 


Ah, 


Ah 


w.s  _ 


C (T  -T  ) 
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sub 
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(68) 


Internal  Heat  Conduction.  X^  is  the  ratio  of  the 
internal  heat  flux  to  the  heat  absorption  rate  due  to 
sublimation.  This  term  is  composed. of  a transient  rate  as 
the  model  cools  and  a steady  state  rate  if  any  temperature 
differences  persist.  The  transient  heat  flux  can  be 
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estimated  with  standard  techniques  and  will  "be  discussed  in 
the  next  chapter.  A mechanism  will  be  described  later 
whereby  the  front  of  the  model  could  be  maintained  at  a 
different  temperature  than  the  sides.  Assuming  that  the 
front  of  the  model  responds  to  this  in  a way  similar  to  a 
sphere  with  a cosine  temperature  distribution  on  one 
hemisphere  and  a ’constant  temperature  on  the  other  (Case 
II),  the  results  of  the  glass  bead  analysis  can  be  applied. 
The  predicted  heat  flux  at  the  high  temperature  point  is 


4r 

kAT 


0.8 


(69) 


Therefore 


qint  _ (0.8)k(^s~  Tint^ 


m Ah 


sub 
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(70) 


sub 


The  Conservation  of  Specie  Equation 

The  result  of  the  more  thorough  analysis  showed  that 
Equation  (17)  of  Chapter  II  also  applies  to  the  ablation 
of  the  composite  material.  In  this  case,  the  concentration 
of  COg  at  the  fluid  dynamic  wall  is  not  necessarily  that 

determined  by  phase  equilibrium  at  the  wall  temperature. 

* • 

Three  phenomena  could  cause  the  concentration  to  be 
otherwise.  The  first  is  the  possibility  that  the  sublima- 
tion temperature  is  different  from  the  wall  temperature, 
the  second  is  that  the  concentration  could  be  attenuated 
between  the  subliming  surface  and  the  fluid  dynamic  wall, 
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•and  finally,  because  of  the  mass  transfer,  the  vapor 
pressure  of  the  C02  at  the  subliming  surface  could  be 
something  less  than  the  equilibrium  value. 

Temperature  Difference . The  possibility  of  a differ- 
ence between  the  fluid  dynamic  wall  temperature,  T , and 

w 

sublimation  temperature,  T , must  be  accounted  for  in  the 
simultaneous  solution  of  the  mass  and  energy  conservation 
equations.  It  is  the  fluid  dynamic  wail  temperature  which 
is  required  to  determine  energy  transfer;  the  sublimation 
temperature  controls  the  mass  transfer.  In  the  simultane- 
ous solution,  the  temperature  difference  is  assumed  to  be 
the  difference  across  the  glass  beads,  and  is  computed 
using  Equation  (65). 

Concentration  Attenuation.  If  the  columns  of  gas 
between  the  glass  beads  were  essentially  stagnant  with 
respect  to  the  boundary  layer  flow,  the  only  mass  flux 
would  be  normal  to  the  surface.  This  can  be  most  reason- 
ably imagined  at  or  very  near  the  stagnation  point  where 
velocities  parallel  to  the  wall  are  quite  low.  The  result 
of  this  situation  is  a diffusion  well  effect  which  causes 

attenuation  of  concentration  levels  (Ref  4s 452).  Given  the 

• * 

* • « 
arbitrary  well  geometry  shown  in  Figure  20,  the  quasi-one 

dimensional  equation  for  the  flow  rate  of  a specie  is 

dK- 

A{x)  [Ki(pv)  - PD  = Mi  (71) 
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If  the  area  is  constant,  the  resulting  concentration 
distribution  is 


1-K 


1-K 


2 = «p(gy) 


(72) 


Furthermore,  if  a surface 
were  completely  covered 
with  one  dimensional 
wells  of  infinitesimal 
wall  .thickness,  then  the 
mass  flux  in  the  well 
(pv)  would  .also  be  the 
fluid  dyhamic  wall  flux, 
m.  Given  walls  of  finite 
thickness , the  one  dimen- . 
sional  analysis  could  be 
used  with  an  adjustment 
to  the  mass  flux  by  the 


Figure  20,  Arbitrary  Geometry 

for 

Diffusion  Analysis 


ratio  of  well  area  to  surface  area,  F, 


(pv)  = xa/ p (73) 

• An  improved  approximation  for  the  glass  bead  simulation  is 
to  perform  a quasi-one  dimensional  analysis  using  the 
idealized  gas  column  geometry  shown  in  Figure  18.  The 
concentration  attenuation  for  a well  depth  y,  no  greater 
than  a,  is 
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(74) 


where  (pv)w  = "the  mass  flux  at  the  top  of  the  well 


Noneguilibrium  Sublimation.  The  pure  C02  models 

exhibited  the  traits  of  a system  experiencing  nonequilibrium 

sublimation,  which  are  a mass  loss  rate  lower  than  the 

analytical  prediction  of  an  equilibrium  model  and  a wall 

* 

temperature  higher  than  the  prediction  (Ref  l) . The 

definition  of  the  equilibrium  vapor  pressure  is  that  pres- 

« 

sure  which  exists  in  a closed  space  over  a liquid  or- solid 

% 

with  no  net  mass  transfer  through  the  phase  interface  and 
with  no  other  substances  present.  Since  the  ablation 
environment  is  quite  different,  nonequilibrium  effects  are 
possible.  This  effect  has  been  investigated  by  many 
authors.  The  direct  result  of  nonequilibrium  vaporization 
is  a decrease  in  partial  pressure  below  the  equilibrium 
value.  For  nonequilibrium  vaporization  (Ref  18  and  23) 


(Pe4.-?)s  = f2  ^ 

where:  a = The  coefficient  of  vaporization 

R = The  universal  gas  constant 
M = Molecular  weight 
and  s indicates  the  phase  interface 
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The  coupled  effect  of  the  inclusion  of  nonequilibrium 
vaporization  in  the  heat  and  mass  transfer  equations  for 
ablation  of  pure  C02  in  air  is  shown  in  Figure  21.  The 
effect  is  negligible  at  low  mass  flux  levels,  but  as  the 
ratio  of  mass  flux  to  coefficient  of  vaporization  reaches  a 
threshold  value  of  10  kg/m  *s,  the  effect  becomes  quite 
dramatic . 

Application  of  the  nonequilibrium  vaporization 
equation  to  the  glass  bead  - C02  model  requires  a correc- 
tion to  the  mass  flux  term.  The  glass  beads  block  the 
surface  area  available  for  vaporization;  therefore, 
continuity  of  mass  requires  an  increase  in  the  mass  flux  at 
the  phase  boundary  which  aggravates  the  nonequilibrium 
condition.  Again,  using  the  glass  bead  surface  model  shown 
in  Figure  18,  the  C02  mass  flux  at  the  subliming  surface  • 
required  to  sustain  the  mass  flux  m at  the  fluid  dynamic 
wall  is 


ms  = 


m 


2 /3 

2^-2  (y/a)  + (y/a); 


(76) 


Substituting  m into  the  nonequilibrium  equation,  the  model 

S 

for  the  glass  bead  C02  surface  is 


2 /6  RT  M, 

s 1 

2/3-2  (y/a)  + (y/a)2 


(77) 
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10  100  1000 
ACTUAL  MASS  FLUX  / COEFFICIENT  OF  VAPORIZATION  (k«/m2**) 

. Effect  of  Nonequilibrium  Sublimation  on  Mass  F^ux  and  Wall  Temperature 


Computation  Scheme 

The  improved  B*  equations  must  be  solved  simultane- 
ously to  determine  the  overall  effect  on  the  C02  mass  loss. 
A B*  prediction  program  was  written  to  include  the 
discussed  effects  and  is  shown  schematically  in  Figure  22. 
The  system  is  no  longer  independent  of  the  heat  transfer 
coefficient,  so  P’eu*e^  must  be  supplied  as  input.  The 
inclusion  of  nonequilibrium  sublimation  necessitates  the 
input  of  a first  estimate  for  the  temperature  at  the  sub- 
liming surface. 
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INPUT  Po  Jo,  y 
ASSUME  B*,TS 


Figure  22. 

I 

;■»  ■ 


OUTPUT 


Iterative  Scheme  to  Predict  B'  and  Model 
Temperatures  For  Glass  Bead  - C02  Composite 
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VI.  Comparison  of  Theoretical  and  Experimental  Results 


While  it  would  he  desirable  to  subject  the  models  to  a 

broad  range  of  heat  transfer  conditions,  the  actual  range 

was  quite  narrow.  Excluding  Run  4,  the  average  value  of 

• 

stagnation  point  Pe^eC^  for  both  pure  CC>2  and  composite 
experiments  was  0.227  t 0.021  kg/m  *s.  This  stagnation 
point  mean  value  and  the  BLIMP  distribution  around  the 
hemispherical  nose,  scaled  to  the  mean  stagnation  value, 
became  the  base  line  for  the  theoretical  results  which 
follow.  Before  attempting  to  predict  the  composite 
response,*  the  pu*  C02  prediction  model  is  examined.  The 
base  line  environmental  conditions  for  the  pure  C0?  results 
were  P = O.976  atm  and  Ta^r  = 280K,  the  average  for  those 
runs. 

Pure  Carbon  Dioxide  Models 

Two  corrections  to  the  simple  theory  which  could 
affect  the  ablation  performance  of  the  pure  C02  models  are 
nonequilibrium  sublimation  and  radiant  heat  flux.  When 
these  phenomena,  Equations  (50)  and  (75) » were  added  to  the 
simplified  ablation  equations  presented  iri  Chapter  II,  a 
coefficient  of  vaporization  equal  to  0.0005  brought  the 
theoretically  predicted  mass  loss  into  closest  agreement 
with  the  experimental  mass  loss  data  and  improved  the 
temperature  agreement  as  shown  in  Figure  23.  Since  more 
confidence  was  placed  on  the  mass  loss  data,  the 
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TEMPERATURE  (K1  1 RELATIVE  MASS  FLUX 


Experimental  Results  of  Pure  C0£  Models— All  Normalized  by 
Equilibrium  Theory  Stagnation  Point  Mass  Loss  Prediction 


vaporization  coefficient  was  chosen  to  fit  the  mass  loss 
data  rather  than  the  temperature  data,  which  would  have 
resulted  in  a coefficient  of  0.0003.  The  effect  of  the 
nonequilibrium  sublimation  is  to  reduce  the  partial 
pressure  of  the  C02  at  the  surface.  Referring  to  the 
conservation  of  species  equation 


M1  Plw/P 

B'  = W2  1 - (Plw/P) 


(19) 


a reduced  partial  pressure  results  in  a reduced  B’.  In 
order  to  then  satisfy  the  energy  equation 


B1 


^r“^2w 

4>Wi + x> 


(W) 


there  must  be  an  increase  in  wall  temperature.  The  radia- 
tion term  is  only  a 1.5$  correction  and  hardly  affected  the 
solution. 

The  evaporation  coefficients  of  most  materials  are  not 
known.  Scala  (Ref  18)  lists  several  of  the  known  coeffi- 
cients which  vary  from  approximately  one  for  water  ice  to 

for  red  phosphorus.  It  is  also  known  that  the  value 
of  the  coefficient  reduces  with  a decrease  in  temperature 
(Ref  9)*  While  a value  of  0.0005  is  lower  than  most  known 
values,  it  does  not  appear  to  be  out  of  the  question. 

There  are  not  many  other  reported  C09  mass  loss  experiments 
to  examine.  Spalding  and  Christie  (Ref  25)  report  experi- 
mental results  which  were  13  to  20  percent  below  their 
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theoretical  predictions;  however,  the  mass  loss  rates  were 
sufficiently  low  that  nonequilibrium  effects  could  only 
improve  their  agreement  by  less  than  two  percent. 

As  a result  of  these  corrections,  the  predicted  B'  and 
T are  no  longer  constant  around  the  body.  The  variable 
wall  temperature  gives  rise  to  steady  state  heat  transfer 
through  the  nosetip;  therefore,  internal  heat  conduction 
should  have  also  been  listed  as  a potential  correction 
term.  However,  the  slight  temperature  difference  for  the 
pure  CC>2  models  is  too  small  to  significantly  affect  the 
solution. 

Glass  Bead  - Carbon  Dioxide  Models 

Nonequilibrium  Model.  Using  the  iterative  scheme 
outlined  in  Figure  22,  ablation  predictions  were  obtained 
for  the  glass  bead-C02  models  using  the  base  line  heat 
transfer  distributions  and  average  environmental  conditions 
of  P = 0.967  atm  and  = 273K.  Parametric  solutions  in 
terms  of  stagnation  point  diffusion  well  depth  were 
obtained  using  the  previously  deduced  value  of  0.0005  for 
the  coefficient  of  vaporization.  In  the  case  of  the 
composite  models,  the  temperature  variation  around  the  body 
is  significant  and  a consistent  method  is  required  to 
determine  the  internal  heat  flux.  As  discussed  in  Chapter  V, 
the  side  wall  temperature,  assumed  equal  to  the  temperature 
of  the  subliming  surface  at  the  90  degree  location,  is  used 
in  conjunction  with  the  Case  II  spherical  conduction 
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solution  to  approximate  the  internal  conduction  term. 
Equation  (70).  The  special  input  conditions  for  the 
solution  at  the  90  degree  location  were  the  convective  heat 
» flux  from  the  BLIMP  distribution,  negligible  internal  heat 
transfer  as  predicted  by  the  Case  II  solution,  and  a 
diffusion  well  depth  which  could  be  different  from  that 
specified  for  the  stagnation  point.  The  average  geometrical 
well  depth  is  probably  the  same  around  the  body  and  was 
assumed  to  be  so  for  the  calculation  of  mass  flux  at  the 
subliming  surface  in  Equation  (77).  Yet,  the  depth  for 
diffusion  purposes  probably  decreases  as  increased  boundary 
layer  velocities  promote  circulation  between  the  beads.  In 
order  to  bracket  this  effect,  the  following  distributions 
were  examined} 

1)  y&iff  = y (A  constant  around  the  hemisphere) 

2)  ydiff  = y (cos  0) 

3)  Ydiff  ” y ) (See  Figure  3) 

e max 

For  each  case,, a geometrical  well  depth,  y,  was  assumed, 

the  ablation  solution  for  the  90  degree  location  was 

obtained,  and  Tg  for  that  location  was  specified  as  the 

. * • 

internal  temperature  for  the  ablation  solutions  around  the 
body. 

The  stagnation  point  solutions' are  shown  in  the  B* 
versus  T plot,  Figure  2k,  for  Cases  2 or  3 which  had  a 
, negligible  well  depth  at  the  90  degree  location.  If  the 
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Figure  24.  Effect  of  Well  Depth  on  Stagnation  Point  Ablation  Response 
of  Composite  Models  with  Nonequi librium  Sublimation  Model. 
( T = 273K,  P = 0.967  atm,  <*  = 0.0005,  Side  Wall  Wells 
Have  Negligible  Diffusion  Depth) 


entire  stagnation  point  mass  flux  decrease  is  attributed  to 
a decrease  in  B',  the  required  B’  value  is  0.126,  which  is 
obtained  with  an  average  geometrical  well  depth  of  50  ixm.  . 
The  effect  of  this  depth  on  mass  loss  and  temperature 
around  the  nose  is  shown  in  Figure  25*  This  model  is  quite 
insensitive  to  the  way  in  which  diffusion  well  depth  is 
varied  about  the  body.  The  mass  loss  distribution  is 
reasonably  well  described  by  these  models;  however,  the 
predicted  internal  temperature  of  186K  for  the  zero  side 
wall  cases  is  higher  than  the  observed  average  asymptotic 
value  of  182K. 

i 

Equilibrium  Model.  Since  the  deduced  value  of  the 
vaporization  coefficient  could  not  be  substantiated  in  the 
literature,  an  equilibrium  sublimation  model  was  also 
investigated.  The  nonequilibrium  effect  can  be  easily 
nullified  by  specifying  a value  of  one  for  the  coefficient 
of  vaporization.  In  this  case  the  assumption  is  made  that 
the  aerodynamic  prediction  of  heat  and  mass  transfer  for 
this  shape  is  6.5%°  high  for  both  the  pure  C02  and  the 
composite  which  brings  the  base  line  stagnation  point  mass 

p 

flux  down  to  0.212  kg/m  *s.  The  average  experimental  mass 
flux  is  23%  below  the  new  base  line,  which  is  equivalent  to 
a reduction  of  B'  to  0.135* 

The  ablation  response  as  a function  of  well  depth  for 
the  cases  with  negligible  side  wall  diffusion  depth  are 
shown  in  Figure  26.  In  order  to  reduce  B'  to  0.135*  a well 
depth  of  80  urn  is  required,  which  is  in  excess  of  the  bead 
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Figure  26.  F.ffect  of  Well  Depth  on  Stagnation  Point  Ablation  Response 
of  Composite  with  Equilibrium  Sublimation  Kodel  ( T * 273K 
P = 0.967  atm,  Side  Wall  Wells  Have  Negligible  Diffusion 


radius  of  60  pm.  The  quasi-one  dimensional  solution  for 
the  concentration  attenuation  probably  under  predicts  the 
actual  two  dimensional  axisymmetrical  process,  so  an 
"effective  well  depth"  of  80  pm  could  have  physical 
significance . 

It  was  possible  to  reduce  the  B'  to  0.135  lor  the  case 

of  constant  well  depth;  however,  an  effective  depth  of 

almost  three  times  the  bead  radius  was  required.  Referring 

to  Figure  26,  the  success  of  the  zero  side  wall  well  depth 

assumption  in  reducing  B*  was  a strong  function  of  the 

internal  heat  conduction  promotjed  by  the  low  side  wall 

r • 

temperature.  The  parameter  X is  the  measure  of  divergence 
of  the  sblution  from  the  simple  ablation  theory  energy 
equation  which  is  also  plotted?  a positive  X results  in  a 
solution  below  the  line,  negative  X solutions  are  above. 
With  a side  well  as  deep  as  the  stagnation  point  well,  the 
internal  temperature  rises  with  well  depth  at  about  the 
same  rate  as  the  stagnation  point  sublimation  temperature, 
thereby  keeping  the  internal  heat  conduction  conti ibution 
to  X low  as  shown  by  the  proximity  of  wall  tempers ture  to 
the  energy  line  in  Figure  27. 

Another  interesting  feature  of  the  model  is  also 
exhibited  in  Figure  2?;  increasing  well  depth  has  a limit 
in  its  ability  to  decrease  B’.  This  figure  was  generated 
by  increasing  well  depth  up  to  60  pm  using  Equation  (7*0. 
then  adding  a multiplicative  factor  to  the  exponent, 
keeping  the  60  pm  value  constant  in  the  trigonometric 
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Figure  27,  Effect  of  Well  Depth  on  Stagnation  Point  Ablation  Response  of 
Composite  With  Equilibrium  Sublimation  Model  and  Side  Wall 
Well  Diffusion  Depth  Equal  to  Stagnation  Point  Depth 
( T = 273K,  P = 0.967  atm) 
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expression.  Hence,  effective  well  depths  in  excess  of 
60  jim  were  generated  and  indicated  by  the  factor  on  the 
chart.  The  value  x 1.5  can  be  interpreted  as  a well  with 
an  effective  depth  1.5  times  that  of  the  ideal  60  jzm  depth. 
When  the  combination  of  mass  flux  and  well  depth  reaches  a 
sufficiently  high  value,  the  concentration  of  C02  at  the 
bottom  of  the*  well  reaches  100$,  The  stagnation  point 
wells  saturate  before  the  side  wall  wells  since  side  wall 
mass  flux  is  lower.  Increasing  depth  no  longer  affects  the 
stagnation  point  wells,  but  finally  saturates  the  side  wall 
wells  which  reach  the  same  temperature  as  the  stagnation 
point  sublimation  surface-the  carbon  dioxide  surface 
temperature  in  equilibrium  with  carbon  dioxide  vapor  with  a 
partial  pressure  equal  to  the  ambient  pressure. 

The  mass  flux  and  wall  temperature  predicted  for  the 
zero  side  wall  well  depth  cases  are  shown  in  Figure  28.  In 
these  solutions  there  is  a significant  dependence  on  the 
manner  in  which  well  depth  is  varied.  The  predicted  mass 
flux  for  the  case  of  well  depth  as  a function  of  edge 
velocity  actually  increases  with  distance  from  the  stagna- 
tion point  to  about  the  45  degree  position,  which  was 
observed  in  some  of  the  experiments. 

Two  Phase  Flow  Effect  on  Transfer  Coefficients.  The 
previous  results  were  computed  assuming  that  the  boundary 
layer  heat  and  mass  transfer  coefficients  were  unaffected 
by  the  presence  of  the  glass  beads.  Reduced  mass  flux  can 
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also  be  the  result  of  a drop  in  transfer  coefficient  or  a 
combination  of  changes  in  both  B’  and  PeueCjj*  A parametric 
study  was  performed  of  the  combinations  of  Peuec|{»  well 
depth,  and  coefficient  of  vaporization  which  would  yield  a 
mass  flux  of  0,0288  kg/m  *s,  which  is  equal  to  a 28#  reduc- 
tion of  the  simple  theory  mass  flux  prediction  at  the  base 
line  condition.  The  more  probable  zero  side  wall  diffusion 
depth  models  were  used.  The  results  are  shown  in  Figure  29; 
also  shown  are  the  predicted  internal  temperatures  associ- 
ated with  the  parametric  combinations.  It  is  evident  from 
this  figure  that  the  knowledge • of  internal  temperature  and 
either  well  depth  or  coefficient  of'  vaporization  is  suffi- 
cient to  determine  the  combination  of  p u C„  and  B’ 

-r  e e H. 
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responsible  for  the  reduced  ma’ss  flux. 

The  most  consistent  explanation-  of  the  experimental 
results  was  obtained  by  using  the  average  recorded  minimum 
temperature , 182K , as  the  internal  temperature  and  the 
previously  deduced  value  of  coefficient  of  vaporization, 
0.0005.  The  intersection  of  these  values  on  Figure  29 
occurs  near  the  40  jim  well  depth  curve  and  corresponds  to 
a transfer  coefficient  of  0.213  kg/m  *s.  The  average  well 
depth,  which  is  two  thirds  of  the  bead  radius,  seems 
reasonable.  The  slight  reduction  in  transfer  coefficient 
is  not  altogether  unexpected  since  the  beads  must  be 
accelerated  at  the  expense  of  boundary  layer  gas  momentum 
which  thickens  the  boundary  layer.  With  higher  velocities 
and,  consequently,  higher  Reynolds  numbers,  the  effect  on 
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Coefficient  Which  Will  Result  in  the  Stagnation  Point  Mass  Flux  0.0288  kg/nr *s  and  the 
Resulting  Internal  Temperatures  ( T = 273K,  P = 0.967  atm.  Side  Wells  have  Negligible 
Diffusion  Depth) 
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transfer  coefficient  could  be  reversed.  The  disturbances 
caused  by  the  solid  particles  and  the  rough  wall  would  be 
expected  to  cause  boundary  layer  transition  to  turbulent 
flow  at  a lower  than  usual  Reynolds  number,  which  would  be 
accompanied  by  the  increased  heat  transfer  associated  with 
turbulent  flow  {Ref  19s 509).  Prior  to  transition,  the 
• rough  wall  has  little  effect  on  the  laminar  transfer 
coefficients  (Ref  19:580). 

It  should  be  recognized  that  the  single  embedded 
thermocouple  does  not  provide  the  important  temperature 
variations  in  time  and  space  required  to  definitely  specify 

1 t 

the  internal  temperature?  yet  the  fact  that  the  recorded 
temperature  remained  for  a significant  period  at  or  very 
near  the  asymptotic  value  suggests  that  the  indicated 
minimum  temperature  is  a good  estimate  of  the  internal 
temperature  for  the  theoretical  model.  Future  experimenta- 
tion could  revise  the  value  for  the  coefficient  of  vapori- 
zation; although,  as  shown  in  Figure  29,  greater  values 
would  require  deeper  and,  therefore,  more  unreasonable  well 
depths.  In  fact  this  theoretical  model  predicts  that 
values  of  coefficient  of  vaporization  greater  than  about 
. 0,0015  are  not  admissible  with  a 182K.  internal  temperature. 

The  optimum  analytical  model  to  correlate  the  data  of 
this  experiment  is  that  the  glass  beads  caused  a 6$ 
reduction  in  the  stagnation  point  boundary  layer  transfer 
coefficient  and  a 23$  drop  in  B*  due  to  the  presence  of 
beads  at  the  surface.  The  effect  of  this  model  on  mass 
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flux  and  temperature  is  compared  to, the  experimental  data 
in  Figure  30  as  a function  of  the  single  phase  flow 
transfer  coefficient.  This  plot  is'  based  on  the  assumption 
that  the  two  phase  flow  effect  has  caused  a constant  6% 
reduction  in  transfer  coefficient,  which  is  probably  a good 
assumption  oyer  the  narrow  range  of  the  experiment. 

Further  experimentation  is  required  to  improve  the  func- 
tional relationship  with  transfer  coefficient. 

Relative  Magnitudes  of  Correction  Terms 

The  actual  value  of  each  correction  varied  with  the 
particular  case  chosen;  however  certain  trends  were  evident. 
Two  of  the  contributions  to  X,  the  increase  in  enthalpy  of 
the  glass  beads  and  the  evolved  gas  from  Tg  .to  -T  , were 
always  quite  negligible  compared  to  the  heat  absorbed  by 
sublimation.  Typically  they  were  only  one  percent  each  for 
the  base  line  solutions.  The  relative  radiation  flux  term 
was  a negative  two  percent  effect  at  the  base  line  condi- 
tion, and  internal  conduction  at  the  stagnation  point 
varied  from  one  percent  to  nearly  four  percent  for  the 
equilibrium  solutions  which  had  the  lowest  side  wall 
temperatures.  Since  the  radiation  term  has  the  opposite 
sign  from  th**  other  contributions,  the  net  correction  to 
the  energy  equation  was  never  very  significant.  This  is 
graphically  illustrated  in  Figures  24,  26,  and  2?  by  the 
fact  that  the  predicted  wall  temperature  is  quite  close  to 
the  reference  line  plot  of  the  simple  energy  equation. 
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TRANSFER  COEFFICIENT  (kg/m2-s) 
figure  30.  Comparison  of  Experimental  Data  and  Theoretically  Predicted 

Stagnation  Point  Mass  Flux  and  Model  Temperature  as  a 

Function  of  Heat  Transfer  Coefficient  ( T = 273K, 

P » 0.967  atm,  cC-  0.0005,  y = 40 ^/m,  Side  Wall  Wells  Have 

Negligible  Diffusion  Depth) 


The  "transient  conduction  was  estimated  using  the 
solutions  presented  in  Schneider  (Ref  20s 2 46).  While  heat 
and  mass  transfer  processes  actually,  beg.in  when  the  model 
is  removed  from  the  mold,  an  upper  bound  for  the  transient 
conduction  can  be  established  if  it  i's  assumed  that  the 
model  temperature  is  initially  uniform  at  the  measured 
internal  temperature  at  the  time  of  model  insertion 
(Figure  13-A)  and  subjected  to  an  instantaneous  change  in 
surface  temperature  to  the  asymptotic  temperature  (Figure 
13-B) . Since  the  area  of  interest  is  .the  hemispherical 
nosetip,  the  solution  for  a sphere  is  used.  The  thermal 
conductivity  and  diffusivity  of  both  pure  C02  and  the 
composite  are  approximately  equal;  therefore  the  internal 
heat  flux  for  all  models  is  approximately  the  same.  The 
result  of  this  analysis  was  that  the  transient  heat  flux  is 
initially  very  high  and  can  remain  a significant  contribu- 
tion to  the  energy  equation  for  approximately  thirty 
seconds.  The  heat  flux  is  toward  the  surface,  a negative 
contribution  to  X,  and  consequently  causes  a higher  mass 
flux.  The  initial  assumptions  for  this  analysis  must  have 
been  too  conservative  since  no  genera]  trend  of  high 
initial  mass  transfer  rates  was  observed.  The  heat  and 
mass  transfer  generated  in  the  free  convection  portion  of 
the  experiment  must  have  established  the  internal  tempera- 
ture gradient  well  enough  that  the  transient  effects  were 
negligible.  Therefore,  transient  heat  conduction  was  not 
included  in  the  improved  modeling  of  this  experiment. 
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The  most  significant  effects,  were  caused  by  altera- 
tions in  the  application  of  the  simple  conservation  of 
species  equation.  Because  of  the  diffusion  wells,  the 
concentration  of  CO^  was  attenuated  69$  from  the  subliming 
surface  to  the  fluid  dynamic  wall  at  the  stagnation  point 
for  the  nonequilibrium  solution.  The  deeper  wells  required 
for  the  equilibrium  solutions  caused  an'  attenuation  of  more 
than  75$»  The  nonequilibrium  sublimation  effect,  which 
caused  a 0.12  atm  reduction  in  C02  partial  pressure  for  the 
pure  C02  predictions,  was  amplified  to  an  O.73  atm  reduc- 
tion because  of  the  higher  local  mass  flux  due  to  the 
blockage  of  surface  area  by  the  glass  beads.  The  two 
phenomena,  diffusion  well  attenuation  and  nonequilibrium 
sublimation,  have  the  similar  net  effect  in  that  they  cause 
a reduction  of  the  available  C02  concentration  at  the  wall 
for  a given  surface  temperature.  This  reduces  the  potential 
B’  in  accordance  with  Equation  (17)  or  (19),  which  results 
in  a higher  wall  temperature  to  satisfy  the  energy  equation. 
The  difference  in  temperature  from  the  subliming  surface  to 
the  fluid  dynamic  wall  was  only  about  three  degrees.  While 
this  had  a negligible  effect  on  the  energy  equation,  it 
decreased  the  convective  heat  flux  to  the  surface  by  about 

*0 

three  percent,  reducing  the  energy  available  for  sublima- 
tion and,  hence,  the  B’. 

General  Applicability  of  Results 

This  research  is  the  first  investigation  of  thermo- 
mechanical ablation  in  which  a known  amount  of  solid 
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material  was  introduced  into  the  boundary  layer  of  an 
ablating  body.  The  solid  addition  resulted  in  significant 
departures  from  pure  thermochemical  ablation  performance. 

The  causes  of  the  different  response  have  been  identified. 

It  is  the  purpose  of  this  section  to  discuss  those 
phenomena  which  have  applicability  to  the  ablation  perform- 
ance of  actual  reentry  vehicle  materials  or  more  general 
applicability. 

The  commonly  used  assumption  that  the  thermomechani- 
cally  removed  material  departs  the  surface  with  the  same 
enthalpy  as  the  material  coming  up  to  the  ablating  surface 
was  shown  to  be  not  always  true,  but  probably  a good 
approximation.  The  mechanically  removed  mass  fraction  in 
the  simulation  was  very  high,  0.71,  which  would  tend  to 
amplify  this  effect.  Yet,  in  spite  of  the  definite 
temperature  increase  of  the  beads,  when  properly  normalized, 
the  increase  in  solid  enthalpy  contributed  only  a 1 % 
correction  to  the  energy  equation. 

The  possibility  of  two  phase  flow  effects  on  the 
boundary  layer  .transfer  coefficients  has  general  applica- 
bility. For  the  conditions  of  the  experiment,  a 6fo 

reduction  in  transfer  coefficient  was  realized.  This  is 

% • 
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generally  less  than  the  precision  of  current  ablation  test 
data,  or  accuracy  of  prediction  capability,  but  must  be 
considered  as  ablation  prediction  and  testing  techniques 
improve.  The  utilization  of  the  techniques  of  this  experi- 
- ment  to  investigate  two  phase  boundary  layer  phenomena  in 
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general  has  promise.  The  data  reduction  and  analytical 
techniques  required  to  determine  the  effects  on  transfer 
coefficient  have  “been  developed  in  this  research. 

The  concentration  attenuation  through  the  diffusion 
wells;  which  had  a large  influence  on  this  experiment, 
probably  also  affects  the  ablation  performance  of  graphites 
and  composite  materials.  Real  materials  are  not  perfectly 
homogeneous;  as  ablation  progresses  the  surfaces  roughen  &*r 
lower  density  matrix  areas  recede  'faster  than  higher 
density  grains  or  re-inforcements  (Ref.  9 and  17) . Micro- 
cracking increases  indepth  porosity  and  some  indepth 
sublimation  occurs  (Ref  9) . The  glass  bead  simulation 
produced ‘diffusion  wells  which  were  easily  modeled  because 

“f 

of  the  nonsubliming  walls;  yet-  every  ablating  rough  wall, 
especially  in  the  vicinity,  of  the-  stagnation  point,  has 
regions  which  are  not  realizing  the  full  benefit  of  the 
convective  currents  of  the  boundary  layer.  Therefore,  the 
diffusion  well  effect  is  probably  present,  to  some  degree, 
whenever  ablation  occurs. 

Of  general  interest  is  the  result  that  C02  exhibits  a 
rather  strong  nonequilibrium  effect.  Since  nc  reference  to 
this  could  be  found  in  the  literature,  this  should  be 
independently  verified.  Embedding  nonsubliming  particles 
in  a subliming  substance  is  a unique  method  of  amplifying 
nonequilibriura  effects  by  blocking-  the  area  available  for 
mass  transfer.  Unfortunately,  since  the  question  of 
nonequilibrium  sublimation  of  graphite  is  of  current 
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interest,*  there  is  no  known  substance  which  can.  perform 

the  same  function  at  the  tempera-cures  of  ablating  graphite 

(greater  than  JJOOt).  Currently,  graphite  specimens  are 

being  arc  jet  tested  over  a large  range' of  transfer 

coefficient  to  determine  the  departure  from  equilibrium 

predictions.  The  surface  roughness  of  ablating  graphite 

probably  tends  to  alleviate  the  nonequilibrium  effect  since 

roughness  increases  the  surface  area  available  for  mass 
transfer. 


* AFML  Contract  F33615-74-C-5094 


with  the  Acurex  Corporation 


VII.  Conclusion 


The  heat  and  mass  transfer  effects  of  thermomechanical 
ablation  were  investigated  by  simulating  this  phenomena 
with  embedded  glass  bead  spheres  in  solid  carbon  dioxide 
(C02)  ablation  mod'els.  Composite  and  pure  C02  models, 
which  were  approximately  25  mm  diameter  hemispherical 
forebodies  with  cylindrical  afterbodies,  were  ablated  in  a 
low  speed  jet  at  approximately  ambient  pressure  and  0°  C 
conditions.  The  pure  C02  models  performed  generally  in  the 
manner  predicted  by  simple  equilibrium  ablation  theory  and 
laminar  heat  transfer  results;  however,  the  stagnation 
point  mass  flux  was  6.2%  lower  than  that  predicted  and  the 
temperature  was  slightly  higher.  These  anomalies  can  be 
explained  as  a nonequilibrium  sublimation  effect,  with  a 
coefficient  of  vaporization  of  0.0005  yielding  the  best  fit 
of  the  experimental  data.  The  effect  of  the  glass  bead 
addition  was  threefolds  (1)  The  C02  mass  flux  was  reduced 
by  2395  at  the  stagnation  point,  which  was  ZQ%  below  the 
simple  ablation  theory  prediction.  (2)  The  C02  mass  flux 
around  the  hemisphere  did  not  monotonically  decrease  with 
distance  from  the  stagnation  point,  but  remained  constant 
or  increased  slightly  with  distance  from  the  stagnation 
point  to  beyond  the  45  degree  location  before  decreasing. 

(3)  The  internal  temperarare  was  increased  approximately  3K. 


A reduction  in  analytically  predicted  mass  flux  can  "be 
accomplished  by  reducing  either  the  transfer  coefficient, 
peueCH»  or  ‘tlie  potential  for  mass  transfer,  B'.  An  analysis 
was  performed  which  considered  the  possibility  of  a reduced 
transfer  coefficient  and  modeled  those  phenomena  which 
could  change  the  predicted  B ' , which  is  generally 
accompanied  by  a change  in  model  temperature.  The  modeled 
effects  were  nonequilibrium  sublimation,  radiative  heat 
transfer,  internal  heat  conduction,  concentration  attenua- 
tion through  the  diffusion  wells  caused  by  the  protrusion 
of  the  densely  packed  glass  beads  through  the  CO^  surface, 
and  the  consequences  of  heat  flux  through  the  glass  beads 
to  the  subliming  surface.  The  average  diffusion  well  depth 
was  left 'as  a free  parameter  in  the  analysis. 

The  most  consistent  explanation  of  the  data  was  that 
in  the  ablation  of  the  composite  models  there  was  a (rf> 
reduction  in  the  boundary  layer  transfer  coefficient  due  to 
two-phase  flow  effects  and  a 2J>%  reduction  in  the  B*.  C.’he 
reduction  in  B’  was  primarily  the  result  of  reduced 
available  C02  at  the  fluid  dynamic  wall  due  to  the  combined 
effects  of  (l)  concentration  attenuation  through  diffusion 
wells  with  an  average  depth  of  2/3  of  the  radius  of  the 
glass  beads,  and  {2)  reduced  surface  area  available  for 
sublima+ion  which  amplified  the  nonequilibrium  sublimation 
effect.  Some  reduction  in  masr,.  flux  was  due  to  the  insula- 
ting effect  of  the  glass  beads  which  maintained  a 3K  degree 
drop  in  temperature  between  the  fluid  dynamic  wall  and  the 
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subliming  surface,  and  the  predicted,  steady  state  internal 
heat  flux  from  the  stagnation  point  to  the  side  walls. 
Radiative  heat  flux,  transient  internal  conduction,  and 
reduced  heat  flux  to  the  C02  due  to  the-  increase  in  glass 
bead  temperature  at  the  surface  were  relatively  insignifi- 
cant. 

This  research  yielded  aspects  of  th'e  ablation  process 
which  are  applicable  to  reentry  vehicle  materials.  The 
diffusion  well  effect  is  probably  present,  to  some  degree, 
for  any  ablating  material.  Two  phase  flow  effects  on  the* 
boundary  layer  transfer  coefficients  should  be  considered 
whenever  a large  fraction  of  solid  material  enters  the 
boundary  layer.  The  commonly  used  assumption  that  the 
’ solids  depart  the  surface  with  the  same  enthalpy  as  the 
material  coming  up  to  the  surface  has  been  shown  to  be  a 
good  approximation. 
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VIII.  Recommended  Improvements  to  the  Experiment 

Further  experimentation  is  recommended  to  confirm  the 
value'  used  for  the  coefficient  of  vaporization,  improve  the 
estimate  of  internal  temperature,  and  further  investigate 
the  two  phase  flow  boundary  layer  phenomena. 

An  independent  estimate  of  the  coefficient  of  vaporiza- 
tion should  be  made  using  existing  physical  chemistry 
procedures  such  as  the  Knudson  Cell  Experiment.  In  lieu  of 
that,  the  problem  can  be  indirectly  addressed  by  performing 
the  ablation  experiment  with  smaller  beads.  The  result  of 
bead  radius  on  B ' , as  predicted  by  eauilibrium  and  non- 
equilibrium  models,  is  shown  in  Figure  31.  The  equilibrium 
model  predicts  a gradual  trend  toward  the  simple  ablation 
model  prediction  as  bead  radius  and,  hence,  well  depth 
decrease.  The  nonequilibrium  solution  predicts  only  a 
slight  change  in  B'  and  Tw  since  the  percentage  of  surface 
area  blocked  by  the  beads  is  not  a function  of  bead  size; 
therefore,  the  nonequilibrium  effect  remains  constant  even 
though  the  diffusion  well  depth  is  reduced. 

Improved  temperature  measuring  capability  is  strongly 
recommended  in  any  future  experiments.  An  increased  number 
of  quick  response  thermocouples  precisely  positioned  with 
protection  from  the  influence  of  the  center  rod  would 
resolve  questions  of  actual  internal  temperature  and 
internal  temperature  gradients  influencing  either  transient 
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figure  31.  The  Effect  of  Bead  Size  on  Normalized  Mass  Flux  and 
Wall  Temperature 
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or  steady  heat  conduction.  . . 

The  possibility  of  reduced  boundary  layer  transfer 
coefficients  due  to  two-phase  flow  effects  could  be  further 
investigated  using  beads  of  varying  density.  Reduced  bead 
material  density  would  certainly  reduce  the  two  phase  flow 
effect  on  the  boundary  layer;  yet  the  net  effect  on 
predicted  B'  and  internal  temperature  would  be  negligible 
since  this  would  only  affect  the  solution  through  the 
corrections  to  the  energy  equation  which  were  quite  insig- 
nificant. 
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APPENDIX  A 

BLIMP 

The  Boundary  Layer  Integral  Matrix  Procedure  (BLIMP) 

was  written  by  Bartlett  and  Kendall  (Ref  2) . This  code  can 

compute  a boundary  layer  along  an  ablating  body  including 

the  nonsimilar  terms  for  a general  chemical  system  and 

incorporates  approximate  formulations  for  mixture  transport 

properties,  including  unequal  mass  diffusion  and  thermal 

diffusion  coefficients  for  all  species.  A laminar  flow 

solution  was  obtained  up  to  the  hemisphere-cylinder  tangent 

using  the  velocity  distribution  shown  in  Figure  3.  BLIMP 

was  run  in  the  mode  which  uses  a complete  surface  mass  and 

energy  balance  for  wall  boundary  conditions  with  the  "quasi- 

steady"  assumption  that  q^  ^ = m (hg  - hQ) , where  q^.  is 

the  heat  conduction  into  the  interior,  and  h is  the 

o 

enthalpy  associated  with  the  asymptotic  interior  tempera- 
ture. In  order  to  insure  that  the  BLIMP  prediction  would 
be  comparable  to  the  approximate  theory,  the  interior 
temperature  assigned  to  the  C02  for  BLIMP  was  wall  tempera- 
ture predicted  by  the  approximate  theory.  The  ”air"  for 
BLIMP  was  a mixture  of  oxygen  and  nitrogen.  The  input 
environmental  conditions  were:  P = 0.976  atm,  T = 273K, 

velocity  = 30*5  m/s,  and  nose  radius  = 10  ram. 
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- -'BLIMP  solves  the  boundary  layer  equations  directly 
including  the  mass  addition  at  the  wall,  so  the  unblown 
heat  transfer  coefficient  and  blowing  correction  are  not 
normally  determined. ' For  purposes  of  Comparison,  BLIMP 
was  run  a second  time  with  no  mass  addition  at  the  wall 
and  in  the  fixed  wall  temperature  mode  using  the  tempera- 
’ ture  computed  with  mass  addition.  The  comparison  for  the 
stagnation  point  is  shown  in  Table  I in  the  body  of  the 
. dissertation.  The  predicted  heat  transfer  coefficient 
around  the  hemisphere  is  shown  in  Figure  A-l.  The 
predicted  concentration  profiled, across  the  stagnation 
point  boundary  layer  are  shown  in  Figure  A-2. 

Afte'r  the  BLIMP  run,  it  was  discovered  that  the  input 
C02  vapor  pressure-temperature*  equilibrium  correlation  was 
not  very  accurate  in  the  temperature  regime  of  interest 
for  the  experiment.  Since  the  primary  purpose  -of  the 
BLIMP  was  the  verification  of  the  simple  theory  and  the 
prediction  of  the  relative  distribution  of  heat  and  mass 
flux  around  the  hemisphere,  it  was  not  necessary  to 
perform  a computation  with  the  improved  correlation. 
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B*  versus  Wall  Temperature  for  CO 0 

in  Air  with  Experimental  Conditions 

Shown 

Experimental  Asymptotic  Temperature  Values  . . , . .. 

Summary  of  Reduced  Data  for  Stagnation  Point 

Mass  Flux  Prediction . 
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Summary  of  Reduced  Data  for  Stagnation  Point  Mass  Flux  Prediction 
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FOR  POINTS  S THRU  12,  Y • 0.006486  * T ♦ 0.07131 
FOR  POINTS  12  THRU  16,  Y . • 0.00SISS  • T ♦ 01017 

IS  . 33  « B3  7S  aa 

ELRPSED  TIME  <5EC> 

13^' 


IK  123 


REXES  I DN  <Cii>  • "I  RECESSION  <CM> 


RBLRT ! DN  RUN  I I 


II 


I i 

t i’ 
*1  \ 

; I 

1 

? 


I .Ml 


.HB 


JUl 


ICOSSION  ON 

lime:  h 


V tf 


Ft*  POINTS  1 THRU  4,  Y - O.OOS022  • T * O.OOIOS 
FOR  »rN7S  9 THRU  12.  Y»  0.00710J  • T -0.12S4- 
FOR  POINTS  12  THRU  16.  Y • 0.00(141  * T * 0.0604 


IS. 


7S 


ELRP5EI>  TIME  <SEC> 


IBS 


ia 


f 


I.Ml 


■SB  1 


.6B  1 


•HB  1 


.a 


RBLRT I DN  RUN  I I 

RECESSION  ON 

LINE  IB 


if 


cf 


d* 


dM 


d' 


FOR  POINTS  1 THRU  4,  Y • 0.005S47  * T ♦ 0.12IS 
FOR  PP'NTS  12  THRU  16,  Y • 0.00SS47  • T 0.1215 


IS 


a 


MS 


BB 


7S 


SB 


IIS 


ia 


ELRP5ED  TIME  < SEC > 


L35 


RECESSION  <Of> 


Run  5 (also  Identified  as  Run  12} 

Simple  Theory  Stagnation  Point  Mass  Flux:  0.00435  g/cm2. 


Position 


Recession  Rate 


Mass  Flux 


cm/s  (xKT) 


g/cm2-.s  (xlO3) 


Normalized 


Stag  Pt 
15  right 
15  left 
30  right 
30  left 
15  right 
15  left 
>0  right 
>0  left 

4.59  ’ 

4.69 

4.75 

4.77- 

4.96 

4.71 

4.31 

4.49 

- 2.85 

2.91 
2.94 
2.96 
3.07 

2.92 
2.67 
2.79 

0.655 
• 0.669 

0.676' 
0.680 
0.705 
0.671 
0.6K 
0.641 

r5  right 
'5  left 

4.94 

3.06 

0.703 

10  right 

• 

10  left 

RBLRT1DN  RUN  12 


RBLRTIDN  RUN  12 


MECESSttN  tM 
LINE  I 


2fl.  lO  SO  00  IBB  12k?  I HO  168 

ELRFGED  TIME  < SEC > 

144 


Run  6 (Also  Identified  as  Run  16) 

Simple  Theory  Stagnation  Point  Mass  Flux: 

Position  Recession  Rate 

— an/s  (xlO3) 


0.00411  g/cm2*s 

Mass  Flux 
g/cm2-s  fxlO3) 


Normal i2ed 


RECESSION  <CM>  * T RECESSION  <CM> 


RECESSION  <CM>  | RECESSION  <CN> 


4 


r 

i 


•SB 


RBLHTION  RUN  56 

REOESEIDM  DN  ■ ■ 

unc  g 


d* 


FOR  POINTS  10  THRU  14.  Y » O.OOJ177  * T ♦ 0.10J2 


ELRP5ED  TIME  <5EC> 


% 


8 


RECESSION  <CM> 


RECESSION  < CM > . T RECESSION  <CM> 


RBLRTION  RUN  IS 


icaaitM  m 

LINES 


70*  MINTS  1‘TIIIU  s,  r • 0.003539  * T - O.OOH 
TO*  MINTS  S THRU  12.  Y . 0.004147  • T - 0.03212 
70*  MINTS  12  THRO  IS.  Y . 0.004401  • T - 0.056S 


2B  HI 

* 

1 HR  60 

m m ins  \m 

» 

• 

ELRPSED  T 1 

i me;**  < sec  > 

••  » 

RBLRT ! ON 

RUN  IS  . 

■ . 

SECESSION  ON 
LIHE  1 1 

. 

. 

JB 

• 

a 

/ 

MR  MINTS  I THRU  3,  Y '•  0.002374  • T ♦ O.OOOS71 
FOR  MINTS  S THRU  10,  Y ■ 0.003444  • " - 0.009S4 
MR  MINTS  11  THRU  16,.  Y • 0.00J163  • t ♦ p.02321 

20  HSI  W 00  • Its  128  | HD 

ELRP5ED  TIME  < SEC  > 


Run  7 (Also  Identified  as  Run  21-1) 

Simple  Theory  Stagnation  Point  Mass  Flux:  0.00392  g/cm  *s 


Position 


Recession  Rate 

_ /„1A 3» 


cm/s  (xlQ- 


Mass  Flux 


Normalized 


Stag  Pt 

4.58 

2.84 

15  right 

4.57 

2.83 

15  left 

4.42 

2.74 

30  right 

4.65 

2.88 

30  left 

- 4.64 

2.88 

45  right 

4.55 

2.82 

45  left 

4.39 

2.72 

60  right 

4.10 

2.54 

60  left 

4.47 

2.77 

75  right 

2.90 

1.80 

75  left 

2.93 

1.82 

90  right  • 

— — 

— — 

90  left 

— 

" i 

0.724 

0.722 

0.699 

0.735 

0.735 

0.719 

0.694 

0.648 

0.707 

0.459 

0.464 


153 


I MB 


ELRP5E1>  TIME*  < SEC  > 
RBLRT 1 DN  RUN  211 

BKE35ICN  ON 

line:  18  . 


RECESSION  «H> 


RBLRT I DN  RUN  2 ! ! 


9CCEH51CN  Ot 
LUCS 


MO  MINTS  1 THOU  7,  Y • 0.002903  * T ♦ 0.000609 
TOO  MINTS  10  THRU  13.  Y • 0.004663  • T - 0.12$ 
MO  MINIS  13  THOU  IS.  Y'-  0.005003  * T - 0.107 


ELRPSED  TIME  <SEC> 


RBLRT 1 ON  RUN  2 I I 


fCOS5IDN  ON 
LIME  11 


TOO  MINTS  1 THOU  7.  Y • 0.002931  • T ♦ 0.000037 
MO  MINTS  10  THRU  17,  Y • 0.004626  * T - 0.126 


eq  « in  a i hd  m 


ELRPSED  TIME  <SEC> 


158 


Run  8 (Also  Identified  as  Run  21-2) 

2 

Simple  Theory  Stagnation  Point  Hass  Flux:  0.00421  g/cm  *s 
Position'  Recession  Rate  Hass  Flux 


Stag  Pt  5.35  .3.32 

15  right  5.21  3.23 

15  left  5.32  3.30 

30  right  5.01  3.10 

30  left  5.05  3.13 

45  right  4.72  2,9% 

45  left  4.75  2.94 

60  right  4.16  2.58 

60  left  4.28  2.66 

•75  right  3.56  2.21 

75  left  3.06  1.89 

90  right  - — 

90  left 


Normalized 


0.789 

0.767 

0.784 

0.736 

0.743 

0.696 

0.698 

0.613 

0.632 

0.525 

0.449 


IS  \ 31  HK 


ts  m 


IS  ISJ 


ELRP5ED  TIME  < SEC > 


RECESSION.  '•  T RECESSION  <<M> 


162 


RECESSION  <CM> 


l 


RBLRTJON  RUN  212 


fCOSEIM  ON 
Lltt  K 


MR  MINTS  1 THRU  11,  Y • 0.0015*0  • T - 0.0050$ 


IK  12 


ELRP5Et>  TIME  <SEC> 


RBLRTIDN  RUN  2J2 


t *B0  4 BECESjlEN  ON 
LINE  II 


MR  MINTS  1 TURD  4,  Y • 0.0030SS  • T * C.001S4 

MR  MINTS  1 THRU  II,  Y » 0.001705  ♦ T • 0.0245 
Y • O.OhSlOJ  • Y - 0*2482 


llllWktlMKl 


ELRP5El>  TIME  <5EC> 


165 


RELRTION  RUN  221 


RECESSION  «N>  ‘ I RECE5SltW  <<M> 


RBLRT I ON  RUN  221 

necxbicn  at 
LIME  s 


km  rollers  i thru  10.  r • 0.002*32  • t ♦ o.ooos** 


■ 

ELRP5ED  TIME  <5EC>  . • 

RBLRT 1 ON 

' . 

RUN  221  !' 

RECESSION  ON 
LINE  II 

• 

• , 

> 

v. 

TO*  P01NT5  1 THRU  10,  Y • 0.003177  • T - 0.00212 

• r 

► ■-«  — t - .,■■■  ■ | ■ ■ KM  -1  ■■  - — » ■ — — ■ ■ ■ > ■ ■ 

ELRP5El>  TIME  <5EC> 


:«xr 


RSLRT10N  RUN  221 

ito^sax  cm 

LIME  S 


FO*  PQIVTS  1 TOKO  10.  Y • 0.0021*0  * T ♦ 0.00*07 
K*  »I«S  1 TH*J  3.  Y • 0.002***  • T - 0.000*0 


IE  » 


Hi  a Bi  71 


ELRPSEt>  TIME  <SEC> 


RBLRTIDN  RUN  221 

wassiQN  a* 

line:  12 


ro*  roiHTS  ] thru  io,  y • o.ooitso  • t - o.oo*** 


X Hi  SR  I N Cfl  7i 


ELRPSED  TINE  < EEC > 


DS/ME/75-1 
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